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ABSTRACT 

We analyze the luminosity function of the globular clusters (GCs) belonging to the early-type galaxies ob- 
served in the ACS Virgo Cluster Survey. We have obtained maximum likelihood estimates for a Gaussian 
representation of the globular cluster luminosity function (GCLF) for 89 galaxies. We have also fit the lumi- 
nosity functions with an "evolved Schechter function", which is meant to reflect the preferential depletion of 
low-mass GCs, primarily by evaporation due to two-body relaxation, from an initial Schechter mass function 
similar to that of young massive clusters in local starbursts and mergers. We find a highly significant trend 
of the GCLF dispersion a with galaxy luminosity, in the sense that the GC systems in smaller galaxies have 
narrower luminosity functions. The GCLF dispersions of our Galaxy and M31 are quantitatively in keeping 
with this trend, and thus the correlation between a and galaxy luminosity would seem more fundamental than 
older notions that the GCLF dispersion depends on Hubble type. We show that this narrowing of the GCLF 
in a Gaussian description is driven by a steepening of the cluster mass function above the classic turnover 
mass, as one moves to lower-luminosity host galaxies. In a Schechter-function description, this is reflected 
by a steady decrease in the value of the exponential cut-off mass scale. We argue that this behavior at the 
high-mass end of the GC mass function is most likely a consequence of systematic variations of the initial 
cluster mass function rather than long-term dynamical evolution. The GCLF turnover mass Mjo is roughly 
constant, at Mjo — (2.2 ±0.4) x 10^ in bright galaxies, but it decreases slightly (by ^35% on average, 
with significant scatter) in dwarf galaxies with Mb.^^x ^ ~18. It could be important to allow for this effect when 
using the GCLF as a distance indicator We show that part, though perhaps not all, of the variation could arise 
from the shorter dynamical friction timescales in less massive galaxies. We probe the variation of the GCLF to 
projected galactocentric radii of 20-35 kpc in the Virgo giants M49 and M87, finding that the turnover point is 
essentially constant over these spatial scales. Our fits of evolved Schechter functions imply average dynamical 
mass losses (A) over a Hubble time that vary more than Mto, and systematically but non-monotonically as a 
function of galaxy luminosity. If the initial GC mass distributions rose steeply towards low masses as we as- 
sume, then these losses fall in the range 2 x IQ^Mq < A < 10^ per GC for all of our galaxies. The trends in 
A are broadly consistent with observed, small variations of the mean GC half-light radius in ACSVCS galax- 
ies, and with rough estimates of the expected scaling of average evaporation rates (galaxy densities) versus 
total luminosity. We agree with previous suggestions that if the full GCLF is to be understood in more detail, 
especially alongside other properties of GC systems, the next generation of GCLF models will have to include 
self-consistent treatments of dynamical evolution inside time-dependent galaxy potentials. 
Subject headings: galaxies: elhptical and lenticular, cD — galaxies: stai" clusters — globular clusters: general 



L INTRODUCTION 

One of the remarkable features of the systems of globu- 
lar clusters (GCs) found around most galaxies is the shape 
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of their luminosity function, or the relative number of GCs 
with any given magnitude. Historically most important has 
been the fact that these distributions always appear to peak, or 
turn over, at a GC absolute magnitude around Mvjo ~ -7.5 
(e.g., Harris 2001), corresponding roughly to a mass ofMjo ^ 
2 X IO^Mq. The near universality of this magnitude/mass 
scale for GCs has motivated the widespread use of the globu- 
lar cluster luminosity function (GCLF) as a distance indicator 
(see Harris 2001; also Ferrarese et al. 2000), and it has also 
posed one of the longest-standing challenges to theories of 
GC formation and evolution. 

In recent years, some significant amount of attention has 
also been paid to the way that GCs are distributed in mass 
around the peak of the GCLF. Traditionally, the full GCLF 
has most often been modeled as a Gaussian distribution in 
magnitude, corresponding to a lognormal distribution of GC 
masses. However, if one focuses only on the distribution of 
GCs above the point where the magnitude distribution turns 
over, it is found that the mass function can usually be de- 
scribed by a power law (Harris & Pudritz 1994), or perhaps 
a Schechter (1976) function (Burkert & Smifli 2000), which 
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is very similar to the mass distributions of giant molecular 
clouds and the young massive star clusters forming in star- 
bursts and galaxy mergers in the local universe (e.g., Zhang & 
Fall 1999). The main difference between ancient GCs and the 
present-day sites of star-cluster formation is then that the mass 
functions of the latter rise steeply upwards towards masses 
much less than Mjo ^ IO^Mq, far exceeding the observed 
frequency of such low-mass GCs. 

There are two main possibilities to explain this fundamen- 
tal difference. The first is that the conditions of star cluster 
formation in the early universe when GCs were assembling 
may have favored the formation of objects with masses in a 
fairly narrow range around ^ (to the exclusion, in 

particular, of much smaller masses). These conditions would 
no longer prevail in the environments forming young clus- 
ters in the nearby universe. Some theoretical models along 
these lines invoke the ^ 10^ Mq Jeans mass at the epoch of 
recombination (Peebles & Dicke 1968), the detailed proper- 
ties of ^ 10^ Mq cold clouds in a two-phase protogalactic 
medium (Fall & Rees 1985), and reionization-driven com- 
pression of the gas in subgalactic (< 10^ Mq) dark-matter ha- 
los (Cen 2001). 

The second possibility is that GCs were in fact born with 
a wide spectrum of masses, like that observed for young star 
clusters, extending from lO^-W Mq down to ~ lO-'-lO'^M© 
or below. A subsequent transformation to the characteristic 
mass function of GCs today could then be effected mainly 
by dynamical processes (relaxation and tidal shocking) that 
are particularly efficient at destroying low-mass clusters over 
the lifetime of a GC system (e.g.. Fall & Rees 1977; Ostriker 
& Gnedin 1997; Fall & Zhang 2001). Some observational 
evidence has been reported for such an evolution in the mass 
functions of young and intermediate-age star clusters (e.g., de 
Grijs, Bastian & Lamers 2003, Goudfrooij et al. 2004). 

If we take the Occam's-razor view that indeed GCs formed 
through substantially the same processes as star clusters to- 
day, then the picture offered by observations of old GCLFs is 
unavoidably one of survivors. There has been some debate as 
to whether it was in fact the long-term dynamical mechanisms 
just mentioned that were mainly responsible for destroying 
large numbers of low-mass globulars, or whether processes 
more related to cluster formation strongly depleted many low- 
mass protoclusters on shorter timescales (Fall & Zhang 2001; 
Vesperini & Zepf 2003). Even the most massive Galactic GCs 
have rather low binding energies £/, < 10^^ erg (McLaughlin 
2000), so that if conditions were not just right, very many 
protoglobular clusters could have been easily destroyed in the 
earliest ^10^ yr of their evolution, through the catastrophic 
mass loss induced by massive-star winds and supernova ex- 
plosions (see, e.g., Kroupa & Boily 2002; Fall, Chandar & 
Whitmore 2005). Furthermore, any clusters that survive this 
earliest mass-loss phase intact but with too low a concentra- 
tion could potentially still dissolve within a relatively short 
time of - lO'^-lO'' yr (Chernoff & Weinberg 1990). Homoge- 
neous observations of large samples of old GCLFs can help 
clarify the relative importance of such early evolution versus 
longer-term dynamical mass loss in the lives of star clusters 
generally. 

The largest previous studies of GCLFs in early-type galax- 
ies were performed with archival HSTAVFPC2 data. Kundu 
& Whitmore (2001a, b) studied the GCLF for 28 elliptical 
and 29 SO galaxies. They concluded that the turnover magni- 
tude of the GCLF is an excellent distance indicator, and that 
the difference in the turnover luminosity between the V and 



/ bands increases with the mean metallicity of the GCs es- 
sentially as expected if the GC systems in most galaxies have 
similar age and mass distributions. Larsen et al. (2001) stud- 
ied the GCLF for 17 nearby early-type galaxies. They fitted 
Student's t distributions separately to the subpopulations of 
metal-rich and metal-poor GCs in each galaxy, and found that 
any difference in the derived turnovers was consistent with 
these subpopulations having similar mass and age distribu- 
tions and the same GCLF turnover mass scale. Larsen et al. 
also fitted power laws to the mass distributions of GCs in the 
range M ~ IO^-IO^'Mq and found they were well described 
by power-law exponents similar to those that fit the mass func- 
tions of young cluster systems. 

In this paper, we study the GCLFs of 89 early-type galax- 
ies observed by HST as part of the ACS Virgo Cluster Sur- 
vey (Cote et al. 2004). This represents the most comprehen- 
sive and homogeneous study of its kind to date. Some of the 
results in this paper are also presented in a companion pa- 
per (Jordan et al. 2006). In the next section, we briefly de- 
scribe our data and present our observed GCLFs in a machine- 
readable table available for download from the electronic edi- 
tion of the Astrophysical Journal. In ^we discuss two dif- 
ferent models that we fit to the GCLFs, and in § |4] we de- 
scribe our (maximum-likelihood) fitting methodology. Sec- 
tion|5]presents the fits themselves, while ^discusses a num- 
ber of trends for various GCLF parameters as a function of 
host galaxy luminosity and touches briefly on the issue of 
GCLF variations within galaxies. In ^we discuss some as- 
pects of our results in the light of ideas about GC formation 
and dynamical evolution, focusing in particular on the rela- 
tion between our data and a model of evaporation-dominated 
GCLF evolution. In ^we conclude. 

2. DATA 

A sample of 100 early-type galaxies in the Virgo cluster 
was observed for the ACS Virgo Cluster Survey (ACSVCS; 
Cote et al. 2004). Each galaxy was imaged in the F475W (~ 
Sloan g) and F850LP (~ Sloan z) bandpasses for a total of 
750 s and 1210 s respectively, and reductions were performed 
as described in Jordan et al. (2004a). These data have been 
used previously to analyze the surface-brightness profiles of 
the galaxies and their nuclei (Ferrarese et al. 2006ab, Cote 
et al. 2006), their surface brightness fluctuations (Mei et al. 
2005ab; 2007), and the properties of their populations of star 
clusters, mainly GCs (Jordan et al. 2004b, Jordan et al. 2005, 
Peng et al. 2006a) but also dwarf-globular transition objects 
(or UCDs, Ha§egan et al. 2005) and diffuse star clusters (Peng 
et al. 2006b). 

One of the main scientific objectives of the ACSVCS is the 
study of the GC systems of the sample galaxies. We have de- 
veloped a procedure by which we select GC candidates from 
the totality of observed sources around each galaxy, discard- 
ing the inevitable foreground stars and background galaxies 
that are contaminants for our purposes. This GC selection 
uses a statistical clustering method, described in detail in an- 
other paper in this series (Jordan et al. 2007, in preparation), 
in which each source in the field of view of each galaxy is 
assigned a probability pcc that it is a GC. Our samples of GC 
candidates are then constructed by selecting all sources that 
have Pec > 0.5. The results of our classification method are 
illustrated in Figure 1 of Peng et al. (2006a). For every GC 
candidate we record the background surface brightness 4 of 
the host galaxy at the position of the candidate, and we mea- 
sure z- and g-band magnitudes and a half-light radius by 
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fitting PSF-convolved King (1966) models to the local light 
distribution of the cluster (Jordan et al. 2005). Photomet- 
ric zeropoints are taken from Sirianni et al. 2005 (see also 
Jordan et al. 2004a), and aperture corrections are applied as 
described by Jordan et al. (2007, in preparation). 

Note that, as part of the ACSVCS we have measured the 
distances to most of our target galaxies using the method 
of surface brightness fluctuations (SBF; Tonry & Schneider 
1988). The reduction procedures for SBF measurements, fea- 
sibility simulations for our observing configuration, and cal- 
ibration have been presented in Mei et al. (2005ab) and the 
distance catalog is presented in another paper in this series 
(Mei et al. 2007). We use these distances in this paper" to 
transform observed GC magnitudes into absolute ones on a 
per galaxy basis whenever we wish to assess GCLF properties 
in physical (i.e., mass-based) terms or need to compare the 
GCLFs of two or more galaxies. While some galaxies have 
larger distances, the average distance modulus that we employ 
is (m-M)o = 31.09±0.03(random)±0.15(systematic), corre- 
sponding to D = 16.5 ±0.1(random)± 1.1 (systematic) Mpc. 

2.1. GCLF Histograms 

There are three main ingredients we need to construct a 
GCLF for any galaxy. First, we have sets of magnitudes, in 
both the z and g bands, for all GC candidates. As mentioned 
above, we generally isolate GC candidates from a list of all 
detected objects by requiring that pcc > 0.5. Note that here 
and throughout, we use g as shorthand to refer to the F475W 
filter, and z denotes F850LP. Also, all GC magnitudes in this 
paper have already been de-reddened (see § 2.7 in Jordan et al. 
2004a for details). 

Second, we have the (in)completeness functions in both 
bandpasses. Our candidate GCs are marginally resolved with 
the ACS, and thus these completeness functions depend not 
only on the GC apparent magnitude m and its position in its 
parent galaxy (through the local background surface bright- 
ness 4), but also on the GC projected half-light radius Rh. 
Separate z- and g-band completeness functions f{m,Ri,,Ii,) < 
1 have therefore been calculated from simulations in which 
we first added simulated GCs with sizes = 1,3,6,10 pc 
and King ( 1 966) concentration parameter c = 1 .5 to actual im- 
ages from the ACSVCS (making sure to avoid sources already 
present), and then reduced the simulated images in an identi- 
cal fashion to the survey data. We next found the fraction of 
artificial sources that were recovered, as a function of input 
magnitude and half-light radius, in each of ten separate bins 
of background light intensity. The final product is a three- 
dimensional look-up table on which we interpolate to obtain 
/ for any arbitrary values of (m,Ri,,Ih). 

Last, we have the expected density of contaminants as a 
function of magnitude for each galaxy, obtained from analy- 
sis of archival ACS images (unassociated with the Virgo Clus- 
ter Survey) of 17 blank, high-latitude control fields, each ob- 
served with both g and z filters to depths greater than in the 
ACSVCS. We "customized" these data to our survey galaxies 
by performing object detections on every control field as if it 
contained each galaxy in turn. This procedure is described in 
more detail in Peng et al. (2006a, their §2.2). The net result is 
17 separate estimates of the number of foreground and back- 
ground objects, as a function of g and z magnitude, expected 
to contaminate the list of candidate GCs in every ACSVCS 

" We use the distances obtained using the polynomial calibration pre- 
sented in Mei et al. 2007. 



field. 

Of the 100 galaxies in the ACSVCS, we restrict our analysis 
to those that have more than 5 probable GCs, as estimated by 
subtracting the total number of expected contaminants from 
the full list of GC candidates for each galaxy. We additionally 
eliminated two galaxies for which we could not usefully con- 
strain the GCLF parameters. This results in a final sample of 
89 galaxies. The GCLF data for these are presented in Table 

m 

The first column of Table [T] is the galaxy ID in the Virgo 
Cluster Catalogue (VCC: Binggeli, Sandage & Tammann 
1985; see Table 1 in Cote et al. 2004 for NGC and Messier 
equivalents). Column (2) contains an apparent z-band mag- 
nitude defining the midpoint of a bin with width given in 
column (3). This binwidth was chosen to be 0.4 for all galax- 
ies. Columns (4)-(6) of the table then give the total number 
Nz.tot of observed sources in this bin; the number A^^.cont of con- 
taminants in the bin as estimated from the average of our 17 
control fields; and the average completeness fraction /, in the 
bin — all applying to the candidate-GC sample defined on the 
basis of our GC probability threshold, pcc > 0.5. Columns 
(7)-(ll) repeat this information for the galaxy's GC candi- 
dates identified in the g band. Columns (12)-(21) are the 
corresponding z- and g-band data for an alternate GC sam- 
ple defined strictly by magnitude cuts and an upper limit of 
Rh < 0."064 ~ 5 pc (which will include the large majority of 
real GCs; Jordan et al. 2005), rather than by relying on our 
Pcc probabilities. This provides a way of checking that se- 
lecting GC candidates by pcc does not introduce any subtle 
biases into the GCLFs (see also ^below). 

The data in Table [T]can be converted to distributions of ab- 
solute GC magnitude by applying the individual galaxy dis- 
tances given in Mei et al. (2007). If they are used to fit model 
GCLFs, it should be by comparing the observed A^tot against 
a predicted (f Nmodei + Ncom) as a function of magnitude. 
This is essentially what we will do here, although we employ 
maximum-likelihood techniques rather than using the binned 
data. However, before describing our model-fitting methodol- 
ogy, we pause first to discuss in some detail the models them- 
selves. We work with two different distributions in this paper: 
one completely standard, and one that is meant to elucidate 
the connections between observed GC mass distributions and 
plausible initial conditions and dynamical evolution histories. 

3. TWO GCLF MODELS 

The term "globular cluster luminosity function" is custom- 
arily used to refer to a directly observed histogram of the 
number of GCs per unit magnitude. We follow this stan- 
dard useage here, and in addition whenever we refer simply 
to a "luminosity function," we in fact mean the GCLF, i.e., 
the distribution of magnitudes again. We denote the magni- 
tude in any arbitrary bandpass by a lower-case m, and thus 
the GCLF is essentially the probability distribution function 
dN /dm. It is not equivalent to the distribution of true GC 
luminosities, since of course m = C-2.5 log L for some con- 
stant C, so dN jdL = (dN / dm)\dm / dL\ has a functional form 
different from that of dN /dm. 

In this paper, when we speak of GC masses, we denote them 
by an upper-case M and we almost always make the assump- 
tion that they are related by a multiplicative constant to GC 
luminosities, such that m = C' -2.5 log M, with C' another 
constant including the logarithm of a mass-to-light ratio (gen- 
erally taken to be the same for all GCs in any one system, as is 
the case in the Milky Way; McLaughlin 2000). We refer to the 
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number of GCs per unit mass, dN jdM, as a "mass function" 
or a "mass distribution." In the literature, it is sometimes also 
called a "mass spectrum." Its relation to the GCLF is 



1 df^ 



dM Md log M 



cxlO" 



dm 



(1) 



As we have already mentioned, most observed GCLFs 
show a "peak" or "turnover" at a cluster magnitude that is 
generally rather similar from galaxy to galaxy. One impor- 
tant consequence of equation ([T]i is that any such feature in 
the GCLF does not correspond to a local maximum in the GC 
mass distribution: if the first derivative of dN /dm with re- 
spect to m vanishes at some magnitude m-io, then the deriva- 
tive of dN /dM with respect to M at the corresponding mass 
scale Mjo is strictly negative, i.e., the mass function still rises 
towards GC masses below the point where the GCLF turns 
over. (More specifically, the logarithmic slope of dN /dM 
at the GCLF turnover point Mto is always exactly -1; see 
McLaughlin 1994 for further discussion.) 

3 . L The Standard Model 

The function most commonly taken to describe GCLFs is 
a Gaussian, which is the easiest way to represent the peaked 
appearance of most luminosity functions in terms of number 
of clusters per unit magnitude. It is thus our first choice to fit 
to each of the observed GCLFs in this paper. Denoting the 
mean GC magnitude [i = (m) and the dispersion = ((« - 
/i)') 'Z^, we have the usual 
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In terms of GC masses, M, this standard distribution corre- 
sponds to a mass function dN /dM = (dN / dm)\dm / dM\ or, 
since m = constant -2.5 logM (assuming a single mass-to- 
light ratio for all clusters in a sample). 
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with (Jm = <T„i/2.5. 

As will be evident in what follows, the GCLFs in a large 
sample such as ours show a variety of detail that is unlikely 
to be conveyed in full by a few-parameter family of distribu- 
tions. But it is also clear that a Gaussian captures some of the 
most basic information we are interested in investigating — 
the mean and the standard deviation of the GC magnitudes 
in a galaxy — with a minimal number of parameters. It is also 
the historical function of choice for GCLF fitting, and in many 
cases the fit is indeed remarkably good. 

Nevertheless, the Gaussian does have some practical lim- 
itations. Seeker (1992) showed that the tails of the GCLF 
in the Milky Way and M31 are heavier than a Gaussian al- 
lows, and he argued that a Student's f distribution (with shape 
parameter ~ 5) gives a better match to the data. More im- 
portantly, the observed GCLFs in our Galaxy and in M3 1 are 
asymmetric about their peak magnitude, a fact which has been 
emphasized most recently by Fall & Zhang (2001). This was 
implicit in the work of McLaughlin (1994), who advocated 
using piecewise power laws to fit the number of GCs per unit 
linear luminosity — or piecewise exponentials to describe the 
usual number of GCs per unit magnitude. Baum et al. (1997) 
used an asymmetric hyperbolic function to fit the strong peak 
and asymmetry in the combined Galactic and M3 1 GCLF. 



However, all of these alternative fitting functions still share 
another shortcoming of the Gaussian, which is that there is no 
theoretical underpinning to it. Moreover, with the exception 
of the power laws in McLaughlin (1994), there is no obvious 
connection with the mass distributions of the young massive 
clusters that form in mergers and starbursts in the local uni- 
verse. We therefore introduce an alternative fitting function — 
based on existing, more detailed studies of initial cluster mass 
functions and their long-term dynamical evolution — to ad- 
dress these issues. 

3.2. An Evolved Schechter Function 
3.2.1. Initial GC Mass Function 

Observations of young star clusters indicate that the number 
of clusters per unit mass is well described by a power law — 
dN /dM cx M"'' with /3 « 2 — or alternatively by a Schechter 
(1976) function with an index of about 2 in its power-law part 
and an exponential cut-off above some large mass scale that 
might vary from galaxy to galaxy (e.g., Gieles et al. 2006a). 
Perhaps the best-observed mass distribution for a young clus- 
ter system is that in the Antennae galaxies, NGC 4038/4039 
(Zhang & Fall 1999). In this specific case, a pure power- 
law form suffices to describe the cluster dN /dM as it is cur- 
rently known; but a Schechter function with an appropri- 
ately high cut-off mass also fits perfectly well. Thus, as- 
summg dN/dM cx M-P exp(-M/Mc), we find from the data 
plotted by Zhang & Fall (1999) that /3 = 2.00 ±0.04 and 
log (Mc/M©) = 6.3^03 for their sample of clusters with ages 
2.5-6.3 Myr; and /3 = 1.92±0.14 and log(M/M,) = 5.9![Jg 
for ages 25-160 Myr 

The mass functions of old globular clusters in the Milky 
Way and M3 1 can also be described by power laws with /3 ~ 2 
for clusters more massive than the GCLF peak (McLaughlin 
1994; McLaughlin & Pudritz 1996; Elmegreen & Efremov 
1997). And the GC mass distributions in large ellipticals fol- 
low power laws over restricted high-mass ranges, although 
here the slopes are somewhat shallower (McLaughlin 1994; 
Harris & Pudritz 1994) and there is clear evidence of cur- 
vature in dN/dM (McLaughlin & Pudritz 1996) that is bet- 
ter described by the exponential cut-off at very high cluster 
masses in a Schechter function (e.g., Burkert & Smith 2000). 
Theoretical models for GC formation, which aim expressly to 
explain these high-mass features of GCLFs and relate them 
to the distributions of younger clusters and molecular clouds, 
have been developed by McLaughlin & Pudritz (1996) and 
Elmegreen & Efremov (1997). 

The important difference between the mass functions of old 
GCs and young massive clusters, then, is not the power-law 
or Schechter-function form of the latter per se; it is the fact 
that the frequency of young clusters continues to rise toward 
the low-mass limits of observations, while the numbers of 
GCs fall well below any extrapolated power-law behavior for 
M < 2 X 10^ Mq, i.e., for clusters fainter than the classic peak 
magnitude of the GCLF. We therefore assume a Schechter 
function, 

-^<xM^^exp(-Mo/M,), (4) 
dMo 

as a description of the initial mass distribution of globular 
clusters generally. We emphasize again that the fixed power 
law of My- at low masses is chosen for compatibility with 
current data on systems of young massive clusters. The vari- 
able cut-off mass Mc is required to match the well observed 
curvature present at M > lO^M© in the mass distributions of 
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old GC systems in large galaxies. This feature is certainly 
allowed by the young cluster data, even if it may not be ex- 
plicitly required by them. 

A strong possibility to explain the difference between such 
an initial distribution and the present-day dN/dM is the pref- 
erential destruction of low-mass globular clusters by a variety 
of dynamical processes acting on Gyr timescales (see Fall & 
Zhang 2001; Vesperini 2000, 2001; and references therein). 

3.2.2. Evolution of the Mass Function 

Fall & Zhang (2001) give a particularly clear recent 
overview of the dynamical processes that act to destroy glob- 
ular clusters on Gyr timescales as they orbit in a fixed galac- 
tic potential. The main destruction mechanisms are dynami- 
cal friction; shock-heating caused by passages through galaxy 
bulges and/or disks; and evaporation as a result of two-body 
relaxation. Only the latter two are important to the devel- 
opment of the low-mass end of the GCLF, since dynamical- 
friction timescales grow rapidly towards low M, as oc M"' . 
(Cluster disruption due to stellar-evolution mass loss does not 
change the shape of the GC mass function if the stellar IMF 
is universal, unless a primordial correlation between cluster 
concentration and mass is invoked; cf. Fall & Zhang 2001 and 
Vesperini &Zepf 2003). 

Tidal shocks drive mass loss on timescales oc phP„, 
where ph (x M/rI is the mean density of a cluster inside its 
half-mass radius, and is the typical time between disk or 
bulge crossings. Evaporation scales rather differently, roughly 

as Tf-yOcM/ . A completely general assessment of the rela- 
tive importance of the two processes can therefore be compli- 
cated. However, tidal shocks are rapidly self-limiting in most 
realistic situations (Gnedin, Lee, & Ostriker 1999): clusters 
with high enough pi, and on orbits that expose them only to 
"slow" and well-separated shocks (i.e., with both the duration 
of individual shocks and the interval Per longer than an in- 

-I /2 

ternal dynamical time, fdyn Ph ) will experience an early, 
sharp increase in ph in response to the first few shocks. There- 
after Tev <C Tsh, and in the long term shock-heating presents a 
second-order correction to the dominant mass loss caused by 
evaporation. Most GCs today, at least in our Galaxy, appear 
to be in this evaporation-driven evolutionary phase (Gnedin et 
al. 1999; see also Figure 1 of Fall & Zhang 2001, and Prieto 
& Gnedin 2006). 

Fall & Zhang (2001) therefore develop a model for the evo- 
lution of the Galactic GCLF that depends largely on evapo- 
ration to erode an initially steep dN/dMo (in fact, they adopt 
a Schechter function, as in eq. H), for one of their fiducial 
cases). They assume — as is fairly standard; e.g., see Ves- 
perini (2000, 2001) — that any cluster roughly conserves its 
mean half-mass density ph as it loses mass, at least after any 
rapid initial adjustments due to stellar mass loss or the first 
tidal shocks, and when the evolution is dominated by evapo- 
ration. The mass-loss rate'^ is then 

^ev = -dM/dt oc M/fi-h oc (M/r\)^I^ ~ constant , (5) 

where fi-h ^ p"' {v^^l^ oc M^I^R^J'^ is the relaxation time at the 
half-mass radius /?/,. Under this assumption, the mass of a GC 
at any age f is just M(f) = Mq - p^yt. For any collection of 

Note that we use /icv to denote the evaporation mass-loss rate in equa- 
tion (5) in order to be consistent with the notation of Fall & Zhang (2001). 
This should not be confused with our use of the symbol fj, — with different 
subscripts — to represent the mean magnitude in a Gaussian description of the 
GCLF. 



clusters with the same density (for example, those on similar 
orbits, if ph is set by tides at a well defined perigalacticon), 
Pf.y is independent of cluster mass as well as time, and if the 
GCs are coeval in addition, then p^^yt is a strict constant. The 
mass function of such a cluster ensemble with any age t is 
then related to the initial one as in equation (11) of Fall 8c 
Zhang (2001): 



dN 



dN 



dM(t) dMo 



dMn 



DM 



dN 
dMn 



(6) 



Thus, simply making the substitution Mo = M + p^yt in the 
functional form of the original GC mass function gives the 
evolved distribution. An initially steep dN/dMo therefore al- 
ways evolves to a flat mass function, dN/dM(t) ^ constant, 
at sufficiently low masses M ^ Pf.yt. As Fall & Zhang show 
for the Milky Way, and as we shall also see for the early-type 
ACSVCS galaxies, this gives a good fit to observed GCLFs if 
the cumulative mass-loss term p^^yt > IO'^Mq by f ~ 13 Gyr 

The key physical element of this argument, as far as the 
GCLF is concerned, is the linear decrease of cluster mass with 
time. While the quickest way to arrive at such a conclusion is 
to follow the logic of Fall & Zhang, as just outlined, there are 
some caveats to be kept in mind. 

Tidal shocks can be much more important than evapora- 
tion for some globulars. In particular, clusters with low den- 
sities and low concentrations (such that shocks significantly 
disturb the cores as well as the halos), and/or those on ec- 
centric orbits with very short intervals between successive 
bulge or disk crossings, may never recover fully from even 
one shock. Rather than re-adjusting quickly to a situation in 
which Tev ^ Tsh, such clusters may be kept out of dynamical 
equilibrium for most of their lives, significantly overflowing 
their nominal tidal radii. Their entire evolution could then be 
strongly shock-driven. This appears to be the case, for exam- 
ple, with the well known Galactic globular, Palomar 5; see 
Dehnen et al. (2004). The extremely low mass and concen- 
tration of Palomar 5 make it highly unusual in comparison to 
the vast majority of known GCs in any galaxy currently, but 
many more clusters like its progenitor may well have existed 
in the past. This then raises the question of whether consider- 
ing evaporation-dominated evolution alone gives a complete 
view of the dynamical re-shaping of the GC mass function. 
Here, however, it is important that the A^-body simulations of 
Dehnen et al. (2004) show that the late-time evolution of even 
the most strongly shock-dominated clusters is still character- 
ized by a closely linear decrease of mass with time: rather 
than ph being conserved in this case, the half-light radius 
is nearly constant in time, and for a given orbit the mass-loss 
rate is M/rsh oc M/p/, oc R^. While the physical reasoning 
changes, the end result for the GCLF of clusters in this phys- 
ical regime is the same as equation ©.'-^ 

The importance of Fall & Zhang's (2001) assumption of 
a constant pi, for evaporation-driven cluster evolution is its 
implication that the mass-loss rate is constant in time. This 
has some direct support from N-body simulations (e.g., Ves- 
perini & Heggie 1997; Baumgardt & Makino 2003). (Note the 
distinction in Baumgardt & Makino between the total cluster 
mass loss and that due only to evaporation; see their Figure 

Another process that may fall in this regime is impulsive shocking due to 
encounters between GCs and massive concentrated objects like giant molec- 
ular clouds; see Lamers & Gieles (2006) for a recent discussion of this. How- 
ever, this is presumably most relevant to clusters orbiting in disks, where the 
shocks can occur in fairly rapid succession. It is not important for the large 
majority of GCs in galaxy halos. 
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6 and related discussion.) But more than this, if evaporation 
is to be primarily responsible for the strong depletion of a 
GC mass function at scales M < ^i^^t, then after a Hubble 
time /iev? must be roughly of order the current GCLF turnover 

mass Mto ~ 2 x 10^ Mq. Since /iev oc p^j^ , this argument ul- 
timately constrains the average density required of the glob- 
ulars, which can of course be checked against data. In addi- 
tion, satisfying observational limits on the (small) variation of 
GCLFs between different subsets of GCs in any one system — 
for example, as a function of galactocentric position — ^puts 
constraints on the allowed distribution of initial (and final) 
GC densities. 

Fall & Zhang calculate in detail the evolution of the Milky 
Way GCLF over a Hubble time under the combined influ- 
ence of stellar evolution (which, as mentioned above, does 
not change the shape of dN /dM except in special circum- 
stances), evaporation, and tidal shocks (which, again, con- 
tribute second-order corrections to the results of evaporation 

1 12 

in their treatment). They relate /iev to GC orbits, by 

assuming that ph is set by tides at the pericenter of a cluster 
orbit in a logarithmic potential with a circular speed of 220 
km s"' . They then find the GC orbital distribution that allows 
both for an average cluster density high enough to give a good 
fit to the GCLF of the Milky Way as a whole, and for a narrow 
enough spread in pi, to reproduce the observed weak variation 
in Mto with Galactocentric radius (e.g., Harris 2001). 

Ultimately, the GC distribution function found by Fall & 
Zhang in this way is too strongly biased towards radial or- 
bits with small pericenters to be compatible with the observed 
kinematics and pi, distributions of globulars in the Milky Way 
and other galaxies — as both they and others (e.g., Vesperini 
et al. 2003) have pointed out. However, Fall & Zhang also 
note that the difficulties at this level of detail do not neces- 
sarily disprove the basic idea that long-term dynamical evo- 
lution is primarily responsible for the present-day shape of 
the GCLF at low masses. The problem may lie instead in 
the specific relation adopted to link the densities, and thus the 
disruption rates, of GCs to their orbital pericenters. In partic- 
ular. Fall & Zhang — along with almost all other studies along 
these lines — assume a spherical and time-independent Galac- 
tic potential. Both assumptions obviously break down in a 
realistic, hierarchical cosmology. Once time-variable galaxy 
potentials are taken properly into account in more sophisti- 
cated simulations, it could still be found that cluster disrup- 
tion on Gyr timescales can both explain the low-mass side of 
GC mass functions and be consistent with related data on the 
present-day cluster orbital properties, ph distributions, and so 
on. Recent work in this vein by Prieto & Gnedin (2006) ap- 
pears promising, though it is not yet decisive. 

We will return to these issues in ^7.11 First, however, we 
describe an analytical form for dN /dM, which combines the 
main idea in Fall & Zhang (2001) — that evaporation causes 
cluster masses to decrease linearly with time — with a plausi- 
ble, Schechter-function form for the initial dN/dMo. We fit 
the evolved function to the GCLF of the Milky Way, to show 
that it provides a good approximation to the fuller, numerical 
models of Fall & Zhang; and then we fit it to our ACSVCS 
data, to produce new empirical constraints for detailed mod- 
eling of the formation and evolution of GC mass functions 
under conditions not specific only to our Galaxy. 

3.2.3. Fitting Functions for dN / dM and the GCLF 



To summarize the discussion above, we assume that the 
mass-loss rate of any globular cluster is constant in time. Fol- 
lowing Fall & Zhang (2001), we expect that this will occur 
naturally if the disruption process most relevant to the GCLF 
in the long term is evaporation, which plausibly conserves the 
average densities p;, of individual clusters inside their half- 
mass radii. Thus, we continue to denote the mass-loss rate 
by /iev. However, it should be recognized that tidal shocks 
can contribute second-order corrections to /iev and may even, 
in some extreme cases, dominate evaporation (though the net 
result arguably could still be a constant total dM/dt). 

For any set of clusters with similar ages t and similar 
Pi, (and on similar orbits, if these significantly affect pi, or 
add tidal-shock contributions to /iev), the cumulative mass 
loss A = /iev? is a constant, so that each cluster has M{t) = 
Mo - A. Combining equation (|4]i for the initial mass distribu- 
tion with equation (|6]l for its evolution then yields an "evolved 
Schechter function" 



dN 



1 



^exp 



M+A 



(7) 



dM iM+Ay-"'''\ M, 

with Mc allowed to vary between galaxies. Once again, A in 
this expression may vary between different sets of GCs, with 
different densities or orbits, in the same galaxy. The detailed 
modeling of Fall & Zhang (2001) takes this explicitly into 
account. But in what follows, we fit equation (|7]i to GC data 
taken from large areas over galaxies, which effectively returns 
an estimate of the average mass loss per cluster over a Hubble 

1 /2 

time. Since /iev oc when evaporation dominates shocks, 
this implicit averaging is essentially done over the distribution 
of GC mean half-mass densities. 

To relate this evolved mass function to the standard obser- 
vational definition of a GCLF — the number of GCs per unit 
magnitude — we write m = C-2.5 logM, S = C-2.5 log A, 
and nic = C-2.5 logM^, where C is related to the solar ab- 
solute magnitude and the typical cluster mass-to-light ratio in 
an appropriate bandpass. The model then reads 



dN 



-^Q—OAim-nic) 



■ exp [-10' 



-0.4(m-m,-)] 



(8) 



dm ^\Q-OA(in-m,)_^ ]^Q-0.4((5-mt)^ 

In both of equations (|7]) and ([8]), the constants of proportion- 
ality required to normalize the distributions can be evaluated 
numerically. 

Figure [T] illustrates the form of the evolved Schecter func- 
tion, in terms of both the mass distribution dN /dM and the 
GCLF dN /dm. (Note that mass M increases to the right along 
the jc-axis in the upper panel, but — as usual — larger M corre- 
sponds to brighter magnitudes m, at the left of the axis in the 
lower panel.) From the equations above, it is clear that the 
mass Mc or the magnitude nic sets the scale of the function, 
while the ratio A/M^. or the magnitude difference {S-nic) 
controls its overall shape. For very small A ^ Mc (faint 
5 ^ mc), the function approaches an unmodified Schechter 
(1976) function. This is drawn in Figure[T]as the bold, broken 
curves that rise unabated toward low cluster masses or faint 
magnitudes. The magnitude mjo at which the GCLF peaks 
in general can be found by setting to zero the derivative of 
equation ^ with respect to m. This yields 

-|^Q-0.8(mTo-«if)_|_ -[Q-0.4(mTo-«if) _|_ jQ-0.4((5-mr)j _ jQ-0.4(i5-mt) _ q 

(9) 

the solution to which corresponds to a mass of 

-(Mc + A) + ^(Mc + A)2 + 4 AM, 
Mto = :r ■ (lU) 
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Fig. 1. — Top: Evolved Schechter mass functions cIN/dM (eq. (7)), for 
various values of the ratio A/Mc, which fixes the shape of the distribution. 
Curves are arbitrarily normalized. The uppermost, broken curve coiresponds 
to A = 0, i.e., a regular Schechter (1976) function with a power-law exponent 
of -2. For non-zero A, t/W/t/M is flat at low masses. Bottom: GChFs dN /dm 
corresponding to the mass functions in the upper panel (see eq. (8]). Curves 
are again arbitrarily normalized, and the parameter controlling the shape is 
the magnitude difference (<5-ffv) = -2.5 log(A/Mc). For any finite {S — mc), 
the GCLF peaks and turns over a t the magnitude mjo given by equation j9) 
(coiTesponding to the mass in eq. 1101 ), and the faint side of the GCLF always 
approaches the limiting shape dN/dm oc lO^ "*'". An'ows mark the turnover 
points of the models shown here. In the limit (S—Mc) — > +oo (i.e. A <C Mc), 
we have that mxo — * 5, while in the limit (S — mc) -oo (large A ^ Mc), 
the turnover mjo mc. For A/Mr > 10 or (S-iric) < -2.5, the GCLF has 
an essentially fixed shape. 



From either of equations ^ or ( fTOl l. or from the sequence 
of curves in Figure [T] it can be seen that when A <C Mc, 
the GCLF peaks at a magnitude nijo ~ S, i.e., the turnover 
reasonably approximates the average cluster mass loss in the 
model (although mxo is formally always fainter than 6). As 
the ratio A/Mf increases, the GCLF turnover initially tracks 
A but eventually approaches an upper limit set by the ex- 
ponential cut-off scale in the mass function: otto — * '«e as 
(S-nic) -oo (A > Mc). 

For any fixed value of A /M^, Figure [T] shows that in the 
limit of low masses, M ^ A, the mass function in equa- 
tion is essentially flat. As Fall & Zhang (2001) first 
pointed out, this is a direct consequence of the assumption 
of a mass-loss rate that is constant in time. It follows gener- 
ically from equation (|6]l above, independently of the specific 
initial GC mass function. At the other extreme, for very high 
masses M ^ A the evolved dN/dM just approaches the as- 
sumed underlying initial function with A = 0. In terms of the 
GCLF, this means that dN/dm tends (always) to an exponen- 
tial, dN/dm oc lO"""*'", at magnitudes much fainter than the 
turnover; and (for initial Schechter function assumed here) to 
the steeper dN/dm cx lO" "*'" exp [-10-"-4 ('«"'«')] for very bright 
magnitudes. The faint half of the GCLF in this model is there- 



fore significantly broader than the bright half. 

Finally, it is worth considering the widths of the GCLFs 
in the lower panel of Figure [T] in more detail. For A = 0, 
the full width at half-maximum (FWHM) of dN/dm is unde- 
fined, since there is no turnover. As the ratio /S./Mc increases 
and a well-defined peak appears in the GCLF, the distribution 
clearly becomes narrower and narrower As we have already 
discussed, even though formally A/M^ can increase with- 
out limit, the turnover magnitude ultimately has a maximum 
brightness mxo ^ nic. Similarly, the FWHM of the GCLF 
approaches a firm lower limit of FWHM ~ 2.66 mag. This 
includes a limiting half width at half-maximum of HWHM ~ 
L59 mag on the faint side of the GCLF, and a smaller 
HWHM ~ 1.07 mag on the bright side. All of these num- 
bers can be obtained from analysis of equation (O by letting 
{5-mc) -oo, i.e., A/M^. +oo. In this Hmit, the GCLF 
approaches a fixed shape and is free only to shift left or right 
depending on the value of ~ »Jto- This limiting shape is al- 
ready essentially achieved with A/M^. = 10 or {5-mc) = -2.5, 
which is plotted in Figure [T] (even though the turnover is still 
about 0.18 mag fainter than in this case). 

As we will see in 35.21 and 36.1.21 the GCLFs observed in 
the ACSVCS are all best fit with A/M^ > 0.1, or (S-mc) < 
2.5 mag. This is the case also in the Milky Way. 

3.3. Comparison with the Milky Way GCLF 

Figure |2] plots the GCLF and the corresponding GC mass 
function in the Milky Way. The upper panel of this figure 
shows the GCLF dN/dm, in terms of clusters per unit abso- 
lute V magnitude, for 143 GCs in the online catalogue of Har- 
ris (1996).'"* (Note again that cluster luminosity and mass in- 
crease to the left in this standard magnitude distribution.) The 
bold, dashed line is the usual Gaussian representation (eq. 121) 
with parameters given by Harris (2001): 

Afv=-7.4±0.1mag; cri/ = 1.15±0.10mag . (11) 

The bold solid curve is our fit of the evolved Schechter func- 
tion in equation (H)), with 

(Sv= -8.0 ±0.3 mag; m^.v = -9.3 ±0.3 mag . (12) 

The lighter, broken line rising steeply towards faint magni- 
tudes is a normal Schechter function with mc as in equation 
(fT2] | but no mass-loss parameter, i.e., (5 — > oo in equation ([8]). 
The shape of this curve is therefore typical of the distribu- 
tion of logarithmic mass for young massive clusters in nearby 
galaxies. 

The lower panel of Figure |2] contains a log-log represen- 
tation of the Galactic GC mass function, dN/dM. To con- 
struct this distribution, we converted the absolute V magni- 
tude of each GC into an equivalent mass by assuming a mass- 
to-light ratio of Ty = 2 Mq Lq for all clusters (as implied 
by population-synthesis models; see McLaughlin & van der 
Marel 2005). The curves here are the mass equivalents of 
those in the upper panel. Thus the bold, dashed curve traces 
equation (O with 

(log(M/Mo)) =5.2 ±0.04; ctm = 0.46±0.04 (13) 

while the solid curve is equation (|7| with 

log(A/M0) = 5.4±O.l ; log(Mc/M0) = 5.9±O.l (14) 

and the lighter broken curve is equation ^ with 
log(Mc/Mo) = 5.9 and A = — again, representative of 
young cluster mass functions. 

''^ http://physwww.mcmasterca/~hanis/mwgc.dat 



8 



JORDAN ET AL. 




absolute V mag 




10= 10" 10=5 1015 

M [M^] 



Fig. 2. — Top: Fits of a Gaussian (dashed curve) and an evolved Schechter 
function (solid curve) to the Milky Way GCLF, expressed as the (normalized) 
number of clusters per unit of absolute V magnitude. The dot-dashed curve 
is a Schechter function with the same value for Mc as the solid curve but with 
mass-loss parameter A set to zero. Bottom: Corresponding observed GC 
mass function clN/dM, and model fits, derived from the GCLF assuming a 
V-band mass-to-light ratio of 2 Mq ^i^'q for all clusters (McLaughlin & van 
der Mai-el 2005). 



Although both model fits to the GCLF are acceptable in 
a statistical sense, the evolved Schechter function yields a 
significantly lower value. This is because of the clear 
asymmetry in the observed GCLF, which appears as a faint- 
ward skew in the top panel and as a failure of the mass 
function dN/dM to decline toward low masses in the bot- 
tom panel. This behavior is described well by the evolved 
Schechter function but is necessarily missed by the Gaussian, 
which systematically underestimates the number of clusters 
withM<3 X lO^M©. 

As a result of this, the best-fit evolved Schechter func- 
tion yields a GCLF peak which is slightly brighter than the 
Gaussian. From the parameters given just above and either 
of equations (|9]l or (fTOl i. we find a turnover magnitude of 
mjo = -7.5 ±0.1 in the V band, some 0.1 mag brighter than 
the Gaussian turnover in equation (fTTl i. The turnover mass 
implied by the evolved Schechter function is thus Mjo ~ 
(1.75 ±0.15) X 10^ M0, just over 10% more massive than the 
Gaussian fit returns. The intrinsically symmetric Gaussian 
model is forced to a fainter or lower-mass turnover in order to 
better fit the relatively stronger low-mass tail of the observed 
GCLF. We find similar offsets in general between the GCLF 
turnovers from the two model fits to our ACSVCS data (see 
g5]and ^below). 

We reiterate that the parameter A in the evolved Schechter 
function represents the average total mass loss per cluster 
(presumably due mostly to evaporation) that is required to 
transform an initial mass function like that of young clusters 



in the local universe, into a typical old GCLF. Both quali- 
tatively and quantitatively, our model fits in Figure |2] corre- 
spond to the various similar plots in Fall & Zhang (2001). In 
fact, the value A ~ (2.5 ±0.5) x 10^ M© obtained here for 
the Milky Way agrees well with the mass losses required by 
Fall & Zhang for their successful models with the second- 
order effects of tidal shocks included. The simple function in 
equation dTji is thus a good approximation to their much fuller 
treatment of the GCLF. 

It is also worth emphasizing just how close A is to the 
GCLF turnover mass scale. This implies that essentially all 
globulars currently found in the faint "half" of the GCLF are 
remnants of substantially larger initial entities. Equivalently, 
any clusters initially less massive than ~2-3 x lO^M© are 
inferred to have disappeared completely from the GC system. 

Despite any difficulties in detail 03.2.2I and Wl-H that 
might remain to be resolved in this evaporation-dominated 
view of the GCLF, and of GC systems in general, it is im- 
portant just to have at hand a fitting formula like the evolved 
Schechter function. In purely phenomenological terms, it fits 
the GCLF of the Milky Way — which is, after all, still the best 
defined over the largest range of cluster masses — at least as 
well as any other function yet tried in the literature. In particu- 
lar, it captures the basic asymmetry of the distribution without 
sacrificing the small number of parameters and the simplic- 
ity of form that have always been the primary strengths of a 
Gaussian description. But at the same time, it is grounded in a 
detailed physical model with well specified input assumptions 
(Fall & Zhang 2001). Fitting it to large datasets, such as that 
afforded by the ACSVCS, thus offers the chance to directly, 
quantitatively, and economically assess the viability of these 
ideas, in much more general terms than has been possible to 
date. 

4. FITTING METHODOLOGY AND TECHNICAL 
CONSIDERATIONS 

4. 1 . Maximum-Likelihood Fitting 

Given either of the models just discussed — or, of course, 
any other — we wish to estimate a set of parameters for the in- 
trinsic GCLF of a cluster sample using the method of max- 
imum likelihood, following an approach similar to that of 
Seeker & Harris (1993). To do so, we make use of all the 
observational material described in ^2.11 

First, we denote the set of GC magnitudes and uncertainties 
in any galaxy, in either the z or the g band, by {m,, e,„.,}. Sec- 
ond, we write the three-dimensional completeness function 
discussed above as f{m,Rit,Ih), which again depends not only 
on GC apparent magnitude but also on a cluster's half-light ra- 
dius and the background ("sky" and galaxy) light intensity at 
the position of the cluster Third, from our 17 control fields we 
are able to estimate the luminosity function of contaminants in 
the field of any ACSVCS galaxy. We call this function b(m), 
and we determine it by constructing a normal-kernel density 
estimate, with bandwidth chosen using cross-validation (see 
Silverman 1986, §§ 2.4, 3.4). Finally, this further allows us 
to estimate the net fractional contamination in the GC sam- 
ple of each galaxy: B = Nc/N, where A^c = (1/17) E,=i ^cj 
with Nc.i the total number of contaminants present in the i-th 
customized control field, and is the total number of all GC 
candidates in the sample. 

Now, given this observational input, we assume that an in- 
trinsic GCLF is described by some function G(m\Q), where 
6 is the set of model parameters to be fitted. The choices 
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for G that we explore in this paper were discussed in detail 
in ^ Thus, for example, for the Gaussian model of equation 
(|2|l, 9 = {/i, CTm}, while for the evolved Schechter function of 
equation (O, 8 = {S^nic}. We further assume that magnitude 
measurement errors are Gaussian distributed, so that — in the 
absence of contamination — the probability of finding an ap- 
parent magnitude m for a GC with given effective radius /?/,, 
galaxy background lb, and magnitude uncertainty e,„ would 
be 

Gj-Cm I e , , 4 , e,„) = A[h(m\ e„) ® G(m \ e)]f(m,Rh,Ih), 

(15) 

where h{m\e„,) = (27re^) '''^exp(-m^/2e^); denotes convo- 
lution; and the normalization A — a function of the GCLF pa- 
rameters 9 and the GC properties R^, lb, and e„, — is fixed 
by requiring that the integral of Gt over the entire magnitude 
range covered by the observations be unity. 

If a fraction B of sources in a galaxy are contaminants, then 
the probability of having a bona fide GC with magnitude m 
(and given/?/,, etc.) is reduced to {1-B)Gt, and thus the like- 
lihood that a set of GCLF model parameters 9 can account for 
A' total objects with observed magnitudes {m,} and properties 

{Rh,i,h,i,^m.i} is 
N 

Cie,B) = II [{l-B)GTimi\e,Rh,i,Ib,i,e„,,d + Bb{md] (16) 
i=i 

in which it is assumed that the luminosity function b{m) of 
contaminants is also normalized. 

For any chosen functional form G(m\Q) of the intrinsic 
GCLF, we specify some initial parameter values 9, compute 
Gt and b for each observed object in a galaxy, and maximize 
on 9 the product in equation (fTST l. In principle, it is possi- 
ble simultaneously to determine the contamination fraction in 
this way, but in practice we found this to be a rather unstable 
procedure (even small inadequacies in the chosen model for 
G can lead to a maximum-likelihood solution that converges 
to quite unreasonable values for B). Thus, we instead made 
direct use of our prior information from the 17 control fields, 
and fixed this fraction to the rather precise average B that we 
have measured for each galaxy. 

The uncertainties in the fitted parameters 9 are estimated 
by using the covariance matrix calculated at the point of max- 
imum likelihood (e.g., Lupton 1993). These uncertainties in- 
clude the effects of possible correlation between the param- 
eters, but they do not include the additional, unavoidable un- 
certainty arising from cosmic variance in the form of b(m) and 
the expected number B of contaminants in any field. As such, 
they constitute lower limits to the total uncertainty. This is not 
a significant issue for GCLF fits to cluster samples combined 
from several galaxies (see below), but it can be important for 
fits to individual galaxies. 

To deal with this, when we fit any individual GC system, we 
re-run our maximum-likelihood algorithm 17 times, each time 
using the background contamination fraction B as estimated 
from a different one of our 17 control fields (versus using B 
from an average of all control fields to obtain the nominal best 
fit). We record the different sets of best-fit GCLF parameters 

In principle there should be another factor multiplying Gt proportional 
to the marginalization over Rj, of the joint GC distribution in m and Ri, times 
an indicator function which is 1 over the area that satisfies poc ^ We 
neglect this factor here, which is justified a posteriori by the agreement of 
results using GC samples constructed using different selection functions. 



obtained in these trials and use the variance in them to evalu- 
ate the additional uncertainty arising from cosmic variance of 
the background contamination. 

4.2. Bias Tests 

Maximum-likelihood estimators are biased in general. It 
is thus important when deriving conclusions to test the bias 
properties of the estimator used, under circumstances similar 
to the ones under study. We have done this specifically for 
the benchmark case of Gaussian fits to the GCLF. After ob- 
taining mean magnitudes and dispersions from our maximum- 
likelihood routine for the 89 ACSVCS galaxies, we analyzed 
20 simulated datasets per galaxy, using the following proce- 
dure. First, we subtracted the number of contaminants BN in 
the galaxy from the total number of GC candidates there, 
to estimate the expected population A^gc of bona fide GCs. 
We then randomly drew a sample of A^gc magnitudes from 
a Gaussian distribution with a mean fi taken to be the fitted 
maximum-likelihood estimate for that galaxy, and a disper- 
sion chosen from a,,, = 0.4,0.7, 1 or 1.3 mag. (We did 5 sim- 
ulations for each of these dispersions, giving the total of 20 
simulations per galaxy.) The randomly generated objects re- 
placed the A^GC objects in that galaxy's sample with the high- 
est Pec values. The values of Ri, and //, of the latter objects 
plus the simulated magnitude are used to determine the com- 
pleteness value / for each source. A uniform random deviate 
is then computed and if that is larger than the value of / the 
source is discarded, a new magnitude drawn from the Gaus- 
sian and the process repeated until the condition is met. In 
this way the effects of completeness are taken into account. 
The maximum-likelihood procedure was finally run on each 
simulated sample and the output parameters compared to the 
input ones. 

The results of these simulations in the z band are summa- 
rized in Figure |3] There may be slight biases in the recovered 
parameters, with (Af7„,/t7„,) w-0.03 and (A/i/cr,,,) «-0.03, 
although there are no significant trends in these average off- 
sets with galaxy luminosity (i.e., sample size). Moreover, the 
statistical significance of these biases is not high (< 3a), and 
so we choose not to correct for them. As a result, it is possi- 
ble that our output best-fit parameters are biased at the level 
of 3% of the GCLF dispersion; but with the possible excep- 
tion only of the most populous GC system (that of M87= VCC 
1316), this turns out always to be smaller than the formal un- 
certainties on the GCLF parameters (see ^5. 11 1. Note that the 
scatter of the retrieved parameters compared with the input 
ones increases towards fainter galaxy magnitudes because the 
candidate-GC sample size is decreasing, and the variance in 
the estimates of both a and /i scales as ~ 1 /N. 

4.3. Effects of Selection Procedure 

As we mentioned in ^ the procedure we used to construct 
a sample of GC candidates for each galaxy involved assigning 
a probability pcc to each source and allowing into the sam- 
ple only those objects with pcc > 0.5. This may influence 
the resulting observed luminosity function and consequently 
affect the derived parameters of any fitted model. In order to 
check that we do not unduly bias our GCLF fitting by this se- 
lection technique, we also constructed alternate candidate-GC 
samples that do not use the selection on pcc but only apply a 
magnitude cut and an u pper limit of Ri, < 5 pc (cf. the second 
half of Table[T]in ^2.1b . The magnitude distributions of such 
samples are free of any selection effects arising from using 
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Fig. 3. — Top: Fractional difference h.a-/a- between input and recovered 
Gaussian dispersion for simulated GCSs with four different Gaussian disper- 
sions assumed (cr = 0.4,0.7, 1.0, 1.3). Bottom: Difference Afi-^/a^ between 
input and recovered Gaussian mean ^ for the same simulated GCSs as in the 
upper panel. 

the pqc values and are useful for testing the robustness of any 
result. Thus, when we fit GCLFs to any of our data, we have 
verified that consistent conclusions are obtained using either 
of our sample definitions. 

4.4. Binned Samples 

While we always perform GCLF fits to individual galaxies, 
some of the fainter systems suffer from small-number statis- 
tics and/or excessive contamination. We thus constructed still 
more GC samples by combining all candidate clusters from as 
many galaxies as required to reach a total sample size above 
some minimum. Going down the list of our target galaxies 
sorted by apparent B magnitude, we accumulate galaxies until 
the expected number of bona fide GCs (i.e., the total number 
of candidates minus the number of contaminants estimated 
from our customized control fields) is > 200. Although many 
of the brighter galaxies satisfy this condition by themselves, 
we refer to the samples defined in this way as "binned" sam- 
ples.'^ There are 24 of them in all, and they are used in ^ 
particularly, to assess trends in GCLF parameters as a function 
of galaxy luminosity without the significant scatter caused by 
the small numbers of GCs in faint systems. 

We excluded 5 galaxies when constructing the binned samples, namely 
VCC 798, VCC 1 192, VCC 1 199, VCC 1297 and VCC 1327. The first was 
excluded due to the presence of a strong excess of diffuse clusters (Peng et al. 
2006b) and the rest because of their proximity to either M87 or M49, making 
their GC systems dominated by those of their giant neighbour; see j|6.3l We 
additionally excluded all galaxies without available SBF distances. 



Our SBF analysis has shown that some of the ACSVCS 
galaxies have distance moduli significantly different from the 
mean (m-M)o = 31.09 mag for Virgo (Mei et al. 2007), and 
thus simply combining the apparent magnitudes of GCs from 
different galaxies with no correction could artificially inflate 
the dispersion of any composite GCLF. To avoid this, we do 
the binning by first using the SBF distances to transform all 
candidate GC luminosities to the value they would have at a 
distance of 31.1 mag (D = 16.5 Mpc). 

5. MODEL FITS 

In this section we present the results of our maximum- 
likelihood fitting of Gaussians and evolved Schechter func- 
tions to the GCLFs in the Virgo Cluster Survey. Recall that 
any alternative model may be fit to the GCLF histograms in 
Table [T] which can be downloaded from the electronic edition 
of the Astrophysical Journal. 

5.1. Gaussian Fits 

The parameter estimates for an intrinsic Gaussian fitted to 
our 89 individual GCLFs are given in Table |2] There we list 
each galaxy's ID number in the VCC and its total apparent 
magnitude Bgai, both taken from Binggeli, Sandage & Tam- 
mann (1985). Following this are the maximum-likelihood 
values of the mean GC magnitude and dispersion and their 
uncertainties in the g-band (pg,(jg), the same quantities in the 

z-band (/ij,,(T,), the fraction B of the sample that is expected 
to be contamination, and the total number of all objects (in- 
cluding contaminants and uncorrected for incompleteness) in 
the galaxy's candidate-GC sample. The last column of Ta- 
ble |2] gives comments on a few galaxies with noteworthy as- 
pects. Note that the uncertainties in the Gaussian parameters 
include contributions from cosmic variance in the shape and 
normalization of the contamination luminosity function b(m) 
(see ED. 

In Figure |4] we present histograms of the observed GCLFs 
along with the best fitting maximum-likelihood models. The 
galaxies are arranged in order of decreasing apparent Bgai 
magnitude (i.e., the same order as in Table |2]i, and there are 
two panels per galaxy: one presenting the z-band data and 
model fits, and one for the g band. The bin width chosen for 
display purposes here is not the same for all galaxies, but fol- 
lows the rule h = 2{IQR)N~^^^, where (IQR) is the interquartile 
range of the magnitude distribution and is the total number 
of objects in each GC sample (Izenman 1991). 

There are four curves drawn in every panel of Figure |4] 
The long-dashed curve is the best-fit intrinsic Gaussian GCLF, 
given by equation (|2]l with the parameters listed in Table |2] 
The dotted curve is this intrinsic model multiplied by the com- 
pleteness function, f(m,Ri,,Ib), after marginalizing the latter 
over the distribution of R/, and 4 for the observed sources 
in each galaxy.'^ The solid gray curve is our kernel-density 
estimate of the expected contaminant luminosity function. Fi- 
nally, the solid black curve is the sum of the solid gray and 
dotted curves; it is the net distribution for which the likeli- 
hood in equation ( fT6l l above is maximized. 

" In order to marginalize f{m,Ri,,Ii,) one needs to know the distributions 
of Ri, and I/, — information which is not available a priori. Using the full 
observed distributions of R;, and 4 is not possible, because they are affected 
by completeness (e.g., faint GCs with large R/, are less likely to be detected). 
We therefore marginaHze / assuming that the underlying distributions in of 
Rl, and 4 are given by the observed distributions for objects satisfying z < 
22.5 and g < 23.7, which gives samples of objects that can be considered 
complete with high confidence, anywhere in any of our galaxies. 
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Fig. 4. — Histograms of the GCLFs for our sample galaxies. For each galaxy we present the z-band and g-band GCLFs side by side. The VCC name and B 
magnitude of the galaxy are indicated in the upper left comer of the left panel, where we also indicate the total number of sources in each histogram and the 
bin-width h used to construct the histogram. Additionally we show the best-fit model (solid black curve), the intrinsic Gaussian component (dashed curve), the 
Gaussian component multiplied by the expected completeness (dotted curve) and a kernel-density estimate of the expected contamination in the sample (solid 
gray curve). The solid black curve is the sum of the solid gray and dotted curves. The galaxies are ordered by decreasing apparent B-band total luminosity, 
reading down from the upper left-hand comer The parameters of the fits are given in Tablem 
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Fig. 5. — (Left) Estimate of Gaussian dispersion in tlie z-band, a-, versus 
tlie same quantity in the g-band, o-g, for tlie GCLFs of our sample galax- 
ies. Uncertainties are Icr. The line marks the one-to-one correspondence 
between these two quantities. (Right) Difference between estimates of Gaus- 
sian means in the g- and z-bands, fig—fi-, versus the mean color (g-z) of the 
GC systems of our sample galaxies. Uncertainties are Icr. The Hne marks the 
one-to-one correspondence between these two quantities. 

Since we have two realizations of the GCLF for every 
galaxy — one in the z band and one in the g band — we are 
able to check the internal consistency of our model parameter 
estimates. Thus, in Figure|5]we compare the measured Gaus- 
sian means and dispersions in the two bands. The left-hand 
panel of this plot shows the scatter of (i, vs. o-g about a line 
of equality, while the right-hand panel shows the difference in 
fitted means (Hg-fi,) vs. the average GC (g-z) color in each 
galaxy (from Peng et al. 2006a), again compared to a line of 
equality. Both cases show excellent agreement between the 
maximum-likelihood results for the two bandpasses. We con- 
clude that the measurements are internally consistent and that 
our uncertainty estimates are reasonable. 

Finally, we also fit Gaussians to our 24 "binned" GC sam- 
ples, constructed by combining the candidates in as many 
galaxies as necessary to reach net sample sizes of at least 
200 (see 34.41 ). The IDs and total magnitudes of the galax- 
ies going into each of these bins are summarized in Table 
[3] along with the best-fit z- and g-band Gaussian parame- 
ters for each binned GCLF and the best-fit parameters for the 
evolved Schechter function discussed in ^(see just below). 
In Figure|6]we display the binned GCLFs in histogram format, 
along with a number of curves representing the maximum- 
likelihood Gaussian fits. The curves in every panel have ex- 
actly the same meaning as in the individual GCLF fits of Fig- 
ure H] We additionally show in this Figure (as the crosses in 
each magnitude bin of each histogram) alternative GCLFs for 
the binned-galaxy samples, obtained by defining GCs on the 
basis of absolute magnitude and an upper limit on the half- 
light radius Rh 

In ^below, we will compare GCLF systematics as a func- 
tion of galaxy properties for these binned samples vs. the fits 
to individual galaxies. We also note here, without showing 
further details, that repeating the exercises of this Section us- 
ing the samples of GC candidates selected only by magnitude 
and Rh, rather than by a pcc criterion, leads to results that are 
consistent in all ways with those we present below. 

5.2. Fits of Evolved Schechter Functions 

We have performed fits of the evolved Schechter function 
in equation ([8]) — or equivalently, the more transparent equa- 
tion (|7| — to the GCLFs of our individual galaxies and binned 
samples. Here we discuss only the results of fitting the 24 

Note that these alternative GCLFs do not have exactly the same numbers 
of objects as the bar histograms corresponding to GC samples defined by 
pac > 0.5. 



binned GC samples, as the results from fitting to all 89 galax- 
ies separately lead to similar conclusions. 

In all these fits, we enforced the constraint that the fitted 
(average) mass loss A be less than ten times the exponential 
cut-off mass scale M^: A/M^ < 10, or {S-nic) > -2.5 in mag- 
nitude terms. This was done because, as was discussed in 33.21 
(see Figure [T]), for such large ratios of A to Mc the evolved 
Schechter function has essentially attained a universal limit- 
ing shape. The likelihood surface then becomes very flat for 
any greater A /M^ and the fitting procedure has difficulty con- 
verging if this parameter is allowed to vary to arbitrarily high 
values. The majority of our evolved Schechter function fits do 
converge to A /Mc values that satisfy our imposed constraint; 
in only one case does the "best-fit" model have the limiting 
A/Mr = 10. 

We show in Figure|2]the binned-sample GCLF histograms, 
along with model curves analogous to those in Figure |6l 
Again, then, the intrinsic evolved-Schechter model GCLF is 
the long-dashed curve; this model multiplied by the marginal- 
ized completeness function is the dotted curve; a kernel- 
density estimate of the contaminant luminosty function is 
shown as the solid gray curve; and the net best-fitting model 
(sum of dotted and solid gray curves) is drawn as a solid black 
curve. Also as in Figure|6] we use crosses in Figure|2]to show 
the GCLFs inferred in every galaxy bin when we define GC 
samples by simple magnitude cuts and /?/, limits, rather than 
by using our pcc probabilities. 

Comparing Figure |7] with Figure |6l it is apparent that an 
evolved Schechter function describes the GCLFs of bright 
galaxies about as well as a Gaussian does. In some of the 
fainter galaxies there is possibly a tendency for the Schechter 
function to overestimate the relative number of faint GCs, but 
it is difficult to assess how serious this might be. The worst 
disagreements between the model fits and the data tend to oc- 
cur in the very faintest extents of the histograms for the hand- 
ful of the faintest galaxy bins at the end of Figure [T] Indeed, 
the largest discrepancies appear at magnitudes where contam- 
inants account for > 50% of the total observed population. 
Any impression of success or failure for any model in these 
extreme regimes of the GCLFs must be tempered by the real- 
ization that the fitting itself is something of a challenge under 
such conditions. 

This is further illustrated by contrasting, in both Figures 
|7] and |6l the GCLFs for cluster samples selected by magni- 
tude and Ri, only (crosses in the figures), to those for sam- 
ples selected on the basis of pcc probabilities (bars). The for- 
mer samples generally tend to put more objects in the faintest 
GCLF bins, an effect particularly apparent in the faintest 
galaxies. The low-mass end of the GCLF for faint galaxies 
is thus not tightly constrained by our observations; there is 
a fundamental uncertainty, due to contamination, that cannot 
be overcome by any selection procedure. (Note that some of 
the more extreme discrepancies between the different GCLF 
definitions — such as in the faintest magnitude bin of BG 20 — 
are due to the presence in some galaxies of a strong excess of 
diffuse clusters that are classified as contaminants when using 
Pcc to construct the sample; see Peng et al. 2006b). But it is 
still worth recalling, in this context, that the "overabundance" 
of low-mass clusters in the evolved Schechter function, vs. a 
Gaussian, is in fact a demonstrably better description of the 
Milky Way GCLF; see Figure |2l 

The fitted magnitude-equivalents S and nic of the mass 
scales A and Mc, in each of the z and g bands, are recorded 
for each of our binned GCLFs in Table [3] In ^we discuss 
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Fig. 6. — Histograms and Gaussian fits to the GCLFs for our binned-galaxy samples. For each sample, named BG« with n = 0,... .23, we present the z-band 
and g-band GCLFs side by side. The identifier of the galaxy bin is indicated in the upper left corner of the left panel, where we also indicate the number N of all 
sources in the histogram (as chosen by requiring poc ^ 0-5) and the bin- width h used when constructing the histograms. In each panel we show the best fitting 
model (solid black curve), the intrinsic Gaussian component (dashed curve), the Gaussian component multiplied by the completeness fraction (dotted curve) and 
a kernel-density estimate of the expected contamination in the sample (soHd gray curve). The solid black curve is the sum of the solid gray and dotted curves. The 
galaxy bins are ordered by decreasing mean apparent B-band luminosity of the galaxies that went into the sample c onstr uction. Crosses in all panels show the 
histograms that result when GC candidates are selected on the basis of cuts in magnitude and half-light radius; see j|4.3l The parameters of all fits are given in 
Tablell] 
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Fig. 7. — Histograms and evolved Schechter function fits to tlie GCLFs for our binned-galaxy samples. For each sample, named BGn with n = 0,... ,23, we 
present the z-band and g-band GCLFs side by side. The identifier of the galaxy bin is indicated in the upper left corner of the left panel, where we also indicate 
the number N of all sources in the histogram (as chosen by requiring poc > 0.5) and the bin-width h used when constructing the histograms. In each panel we 
show the best fitting model (solid black curve), the intrinsic evolved Schechter component (dashed curve), the evolved Schechter component multiplied by the 
completeness fraction (dotted curve) and a kernel-density estimate of the expected contamination in the sample (solid gray curve). The solid black curve is the 
sum of the solid gray and dotted curves. The galaxy bins are ordered by decreasing mean apparent B-band luminosity of the galaxies that went into the sample 
construction. Crosses in all panels show the histograms that result when GC candidates are selected on the basis of cuts in magnitude and half-light radius; see 
j|4.3l The parameters of all fits are given in Table[3] 
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Fig. 8. — Upper panel: Comparison of the absolute z-band magnitude 
of the GCLF turnover, as inferred from the maximum-likelihood fitting of 
intrinsic evolved Schechter functions versus that inferred from Gaussian fits. 
Lower panel: Comparison of the FWHM of the intrinsic z-band GCLFs. as 
returned by the Gaussian and evolved Schechter fits. Similar plots for the fits 
to our g-band GCLFs look the same as these z-band results. 



in detail the conversion of these to masses and also consider 
dependences of A and on galaxy luminosity. 

Just before looking at these issues, Figure |8] compares 
the turnover magnitudes and full-widths at half maximum 
(FWHM) for the binned z-band GCLFs as returned by the fits 
of evolved Schechter functions (see eq. ||9)), against the same 
quantities implied by our Gaussian fits. For the turnovers, 
there is a slight offset, in that the fitted Schechter functions 
tend to peak at slightly brighter magnitudes (typical differ- 
ence <0.1 mag, corresponding to a turnover mass scale that 
is < 10% larger than implied by the Gaussian fits). This is 
very similar to the offset in the two fitted turnover magnitudes 
for the Milky Way GCLF in ^3.31 As we discussed there, the 
discrepancy is a result of the intrinsic symmetry assumed in 
the Gaussian model, vs. the faint-end asymmetry built into the 
evolved Schechter function. 

The FWHMs differ more substantially between the two 
functional forms, with the evolved-Schechter fits being typ- 
ically ~0.5 mag broader (or about 0.2 dex in terms of mass) 
than the Gaussian fits. But this is again only to be expected 
from the asymmetry of the former function vs. the symmetry 
of the Gaussian. As was noted at the end of 33.21 the shape 
of the evolved Schechter function is universally flat in terms 



of dN/dM for low GC masses, or universally cx lO"'' "*™ in 
terms of dN /dm for magnitudes much fainter than the peak 
of the GCLF. As a result, the faint side of the GCLF is always 
broader than any Gaussian, and so if the two models give com- 
parable descriptions of the bright halves of all GCLFs, the 
FWHM of the evolved Schechter functions must always be 
larger than those of the Gaussian fits. Moreover, for very nar- 
row observed GCLFs, fit by small Gaussian (t,„ (primarily to 
reproduce the steepness of the bright side of the GCLF, as dis- 
cussed below), the evolved Schechter function fits are limited 
by a minimum FWHM of ~2.66 mag ( 33.2b . explaining the 
tendency towards a plateau at the left side of the lower panel 
of Figure m 

6. TRENDS BETWEEN AND WITHIN GALAXIES 

Having fitted two different GCLF models to each of our in- 
dividual galaxies and binned samples, we now outline some 
systematic variations in the properties of GC mass distribu- 
tions indicated by this work. First, we examine the depen- 
dence of GCLF parameters on host galaxy luminosity; then — 
even though the ACS VCS data are not ideal for this purpose — 
we look for any evidence of GCLF trends with radius inside 
die two brightest Virgo galaxies, M49 (VCC 1226) and M87 
(VCC 1316). 

6.1. Variations with Galaxy Luminosity 

6.1.1. Gaussian Parameters 

Figure|9]shows one of the main results of this paper: GCLFs 
are narrower in lower-luminosity galaxies (see also Jordan et 
al. 2006). 

The upper panel of this figure plots the Gaussian disper- 
sion that best fits the z-band GCLF, as a function of absolute 
galaxy magnitude M^^gai for our 89 individual galaxies. Filled 
circles represent galaxies with measured (SBF) distance mod- 
uli, while open triangles correspond to galaxies for which 
no distance modulus is available and for which we assume 
(m-M)o = 31.1 (consistent with the average Virgo distance 
modulus of Mei et al. 2007) to compute Mb.shI- The lower 
panel shows the analogous result for our §-band GC data. The 
straight lines drawn in the panels are convenient linear char- 
acterizations of the (t„-Mb.!,m trends: 

cr, = (1.12±0.01)-(0.093±0.006)(MB,gai+20) (17) 

and 

o-g = (1.14±0.01)-(0.100±0.007)(MB,gai+20). (18) 

While it has been reported before that there is a tendency for 
the GCLFs in lower luminosity galaxies to show somewhat 
lower dispersions (e.g., Kundu & Whitmore 2001a), the ho- 
mogeneity of our sample and analysis make this the most con- 
vincing demonstration to date of the existence of a continous 
trend over a factor of w400 in galaxy luminosity. It is partic- 
ularly noteworthy that the fainter galaxies in our sample — all 
of which are early type — ^have very modest a„, < 1, values 
more usually associated with the GCLFs of late-type galax- 
ies. In fact we have also plotted on Figure |9] the V-band 
GCLF dispersions (Harris 2001) and absolute bulge luminosi- 
ties of the Milky Way (large fifled star at MB.gM = -18.8; de 
Vaucouleurs & Pence 1978) and M31 (large filled triangle at 
■^^B.gai = -19.2; from Kent 1989, but assuming a distance of 
810 kpc). Clearly these fall well in the midst of our new data, 
and thus the correlation of a with M^ gai would appear to be 
more fundamental than the older view, that GCLF dispersions 
depend on galaxy Hubble type (Harris 1991). 
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Fig. 9. — Top: GCLF dispersion a., inferred from Gaussian fits to the 
z-band data, versus galaxy Mg.gai- Filled symbols are galaxies for which 
we have available SBF distances while open triangles represent galaxies for 
which we do not and for which we have assumed a distance modulus of (m — 
M)o = 31.1. Th e dashed Une is the linear relation between cr- and Mg 
in equation U7t . Bottom: Same comparison, but for the intri nsic Gaussian 
dispersion of the ^-band GCLFs, erg. Dashed line is equation )18t . In both 
panels the star shows values for the MW and the triangle represents M3 1 . The 
oudiers at Mg.gA - -21.2 and Mb.^oX ~ -19.9 in both panels are VCC 798 
and VCC 2095, galaxies which have an excess of faint, diffuse star clusters 
(Peng et al. 2006b). 

At this point it should be noted that the GCs in brighter 
galaxies are known to have broader color distributions, and 
hence larger dispersions in metallicity, than those in fainter 
galaxies (e.g., Peng et al. 2006a). But cluster mass-to-light 
ratios, T, are functions of [Fe/H] in general, so there will be 
some galaxy-dependent spread in their values. Since the vari- 
ance in an observed luminosity distribution is related to that 
in the mass distribution, by the usual <J^(\og L) = cr^(log M) + 
cr^(log T), this then suggests the possibility that the trend we 
see in the GCLF a- and (jg vs. galaxy luminosity might re- 
sult from systematics in cr(log T) vs. M^ gai on top of a more 
nearly constant CT(log M). In fact, this idea was recently in- 
voked by Waters et al. (2006) as a potential explanation for 
the fact that the /-band GCLF of M87 is broader than that 
of the Milky Way; and by Strader et al. (2006) as a possible 
reason for the narrower composite GCLF of a subsample of 



ACSVCS dwarfs versus the GCLFs of Virgo giants. How- 
ever, neither of those works checked these claims quantita- 
tively. We have done so here (see also Jordan et al. 2006), and 
we find that the explanation is not tenable. 

As we will discuss further in ^6.L2I GC mass-to-light ra- 
tios in the longer- wavelength z band vary by less than ±10% 
over the entire range -2 < [Fe/H] < 0, which includes the 
large majority of clusters. Thus cr(log T,) < 0.04 no mat- 
ter what the details of the GC metallicity distribution are — 
making for an utterly negligible "correction" to the observed 
(^(logL ) = a,/ 2.5 for all of our GCLFs. In the shorter- 
wavelength g band, mass-to-light ratios are more sensitive to 
cluster colors. But here the close agreement of our g- and 
z-band GCLF dispersions shows immediately that the former 
must be reflecting the properties of the GC mass functions just 
as closely as the latter are. Indeed, more detailed calculations, 
which include the observed specifics of the color distributions 
in our galaxies (Peng et al. 2006a), confirm that the spread in 
expected GC values contributes ~ 0.02 mag to the total 
observed GCLF dispersion — an amount well within the ob- 
servational uncertainties on <jg in the first place'^. Thus, we 
proceed knowing that the correlations between GCLF disper- 
sion and galaxy luminosity that we are discussing here are 
very accurate reflections of equivalent trends in the more fun- 
damental GC mass distributions. 

Because of the symmetry assumed in the model, the trend 
of decreasing Gaussian cr,,, in Figure |9] might appear to imply 
a steepening of the GCLF on both sides of the turnover mass. 
However, as we have already discussed, if we take the more 
physically based, evolved-Schechter function of equation (O 
or (|7]) to describe the distribution of GC masses, then all 
GCLFs must have the same basic shape (and thus half-width) 
for clusters fainter than about the turnover magnitude — in 
which case the trends in Figure |9] can only be driven by sys- 
tematics in the bright side of the GCLF. Indeed, as was men- 
tioned in ^above (and discussed at length by, e.g., McLaugh- 
lin & Pudritz 1996), it has long been clear that power-law rep- 
resentations of the GC mass function above the turnover mass 
in the Milky Way and M3 1 are significantly steeper than those 
in M87, M49, and other bright ellipticals; there is no "univer- 
sal" power-law slope for present-day GC mass functions. 

Given these points, we have also performed maximum- 
likelihood fits of pure power-law mass distribu- 
tions (dN/dM cx M~^; or, in terms of magnitude, 
dN/dm oc 10"4(/3-i>") to GCs between ~ 0.5-2.5 mag 
brighter than the turnover magnitude in the cluster samples 
of our individual galaxies. (Such subsamples are both 
highly complete and essentially uncontaminated in all of 
our galaxies). The best-fit (3 for the 66 galaxies in which 
we were able to measure it are presented in Table H) The 
results from fitting to the z- and g-band data are similar, and 
thus we show only the former here, in the upper panel of 
Figure[TO] This confirms that the high-mass end of the GCLF 
steepens systematically for decreasing galaxy luminosity, 
independently of how the low-mass GC distribution behaves. 
In Figure [TO] we also plot a star and triangle showing (3 
values for the Milky Way and M31 respectively, measured 
in the same mass regime using the data from Harris (1996) 
and Reed et al. (1994) assuming a V-band mass-to-light 
ratio M/Ly = 2. The lower panel of Figure [TO] then plots the 
fitted power-law exponent for high GC masses against the 

" We note that the median value of (ffg — o"z) for our sample galaxies is 
0.02 mag. 
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Fig. 10. — Upper panel shows the slope (3- of the power law that best fits 
the z-band GCLF data for GC masses 3 X 10^ < (M/Mq) < 2 x 10*, against 
host galaxy absolute B magnitude. The star and triangle show /3 values for 
the Milky Way and M3 1 respectively, measured in the same mass regime us- 
ing the data from Hanis (1996) and Reed et al. (1994) assuming a V-band 
mass-to-light ratio M/Lv = 2 Mq Ly^Q- Lower panel shows the correlation 
between this power-law index and the dispersion cr, in a Gaussian represen- 
tation of the GCLF. These graphs illustrate that the systematic "nar rowi ng" 
of the GCLF for decreasing galaxy luminosity, as seen in Figs.l9]and ll2l is a 
real phenomenon rather than an artifact of the Gaussian model: it shows up 
clearly as a steepening of the (largely complete and relatively contamination- 
free) high-mass end of observed GCLFs. Corresponding plots for the g-band 
GCLFs are very similar to these. 

Gaussian GCLF dispersion from Figure |9] showing that there 
is indeed a clear correlation between these two parameters 
in the sense that a narrower Gaussian a reflects a steeper 
high-mass power-law (5. 

The regularity and the high significance of the narrowing of 
the GCLF as a function of galaxy luminosity — or the steepen- 
ing of the mass distribution above the classic turnover point — 
places a new and stringent constraint on theories of the for- 
mation and evolution of the mass function of GCs. In one 
sense, this is then on a par with the modest amount of varia- 
tion seen in the turnover mass. An important difference may 
be that the GCLF turnover could be imprinted to some large 
extent by long-term dynamical evolution (Fall & Zhang 2001 ; 
though see, e.g., Vesperini 2000, 2001, and Vesperini & Zepf 
2003 for a differing view, and ^7. 11 below for a discussion 
of caveats). By contrast, most analyses agree that the shape 
of dN /dm above the turnover is largely resistant to change 



by dynamical processes 07.2l i — in which case it seems most 
likely that the systematic variations in Figure s|9] and [10] are re- 
flecting a fundamental tendency to form massive star clusters 
in greater relative numbers in more massive galaxies. 

Moving now to the GCLF turnover magnitude, in Figure [TTI 
we show the absolute ji, and /ig as functions of host galaxy ab- 
solute magnitude Mfi^gai- In both panels of this figure, horizon- 
tal lines are drawn at the levels of the typical turnovers in large 
ellipticals: excluding VCC 798, which has an anomalously 
large excess of faint, diffuse star clusters (Peng et al. 2006b), 
the average Gaussian turnovers for ACSVCS galaxies with 
Mfi.gai < -18 are 

(/i-) = -8.4±0.2 

(£) = -7.2±0.2. (A^i^,gai<-18) (19) 

The turnover in the Milky Way is shown as a large filled star 
and that in M3 1 is represented by a large filled triangle, as in 
Figure |9] We estimated these turnovers from the V-band val- 
ues given in Table 13 of Harris (2001), by applying (g-V) and 
{V-z) colors calculated for 13-Gyr old clusters with [Fe/H] = 
-1.4 for the Milky Way (Harris 2001) and [Fe/H] = -1.2 for 
M31 (Barmby et al. 2000) using the PEGASE population- 
synthesis model (Fioc & Rocca-Volmerange 1997). 

The z-band turnovers in the upper panel of Figure [TTI show 
a tendency to scatter systematically above (fainter than) the 
bright-galaxy value for systems with M^.gai ^ -18, but there 
are no such systematics in the g-band turnovers in the lower 
panel. Interpreting these results is most easily done in terms 
of equivalent turnover mass scales, and thus we defer further 
discussion to ^6.1.21 where we use the PEGASE model to 
convert all of our GCLF parameters to their mass equivalents. 
We note here, however, that the near constancy of Hg in Figure 
[TTl is equivalent to the well known "universality" of the GCLF 
turnover in the more commonly used V band (since our g is 
die HST F475W filter, which is close to standard V). 

Before discussing masses in detail, we plot in Figure[T2]the 
Gaussian means and dispersions of the z-band GCLFs in our 
24 binned samples, vs. the average absolute magnitude of the 
galaxies in each bin (see Table [?]i. The straight lines in each 
panel are just those from the upper panels of Figs. [TTI and [9l 
characterizing the fits to all 89 individual galaxies. This com- 
parison shows that the results from our single- and binned- 
galaxy GC samples are completely consistent, so that our bin- 
ning process has served — as intended — to decrease the scatter 
in the observed behavior of /i and a at low galaxy luminosi- 
ties. It also confirms the results of our simulations in ^4.2! 
above, which showed that our maximum-likelihood model fit- 
ting is not significantly biased by size-of-sample effects. A 
plot like Figure[T2] but using our Gaussian fits to the individ- 
ual and binned ^-band GCLFs, leads to the same conclusions. 

6.1.2. Mass Scales 

To better understand the GCLF trends discussed above, and 
to mesh the Gaussian-based results with those from fits of the 
more physically motivated evolved-Schechter function, it is 
advantageous to work in terms of GC mass, rather than z and 
g magnitudes. To make this switch, we rely on population- 
synthesis model calculations of ig-z) colors and g- and z- 
band mass-to-light ratios as functions of metalUcity for "sim- 
ple" (single-burst) stellar populations. 

The model we use is version 2.0 of the PEGASE code (Fioc 
& Rocca-Volmerange 1997), which we have run by inputting 
the stellar initial mass function of Kennicutt (1983) to com- 
pute cluster masses and g and z luminosities as functions of 
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Fig. 11. — Top: Absolute magnitude fi- of the GCLF turnover, versu.s 
Mg gal, inferred from Gaussian fits to the z-band GCLFs. Filled symbols are 
galaxies for which we have available SBF distances while open triangles rep- 
resent galaxies for which we do not and for which we have assumed a distance 
modulus of {m—M)o = 31.1. The dashed line is at fiz = —8.4, the average for 
galaxies brighter than Mg.a-A = —18. Fainter galaxies have turnover magni- 
tudes that tend to scatter fainter than this. Bottom: Same comparison, for 
the g-band magnitude of the GCLF turnover. Horizontal line is fig = -7.2. 
In both panels the star shows values for the MW and the triangle represents 
M31. The outliers at Mg gai ~ -21.2 and Mfi gai ~ -19.9 in both panels are 
VCC 798 and VCC 2095, galaxies which have an excess of faint, diffuse star 
clusters (Peng et al. 2006b). 

age for several fixed values of [Fe/H]. The results, at an as- 
sumed uniform GC age of 13 Gyr, are illustrated in Figure[T3l 
which plots the mass-to-light ratios and T, in solar units 
and the (g-z) color against [Fe/H]. Given the average (g-z) 
of the GCs in any of our galaxies (from Peng et al. 2006a), we 
interpolate on these PEGASE model curves to estimate aver- 
age g and z mass-to-light ratios. Table |5] lists the mean GC 
color in each galaxy and the M/L values we have derived. 

It is clear from Figure[T3]and Table|5]that the z-band mass- 
to-light ratio varies by only a modest amount for most GCs 
in our samples: we generally have 0.8 < {(g-z)) < 1.2 in 
these cluster systems, and thus 1.45 ^ T- < 1.55 MqLq. 
A z-band luminosity is therefore a very good proxy for to- 
tal cluster mass. By contrast, over the same range of GC color 
or metallicity, the g-band mass-to-Ught ratio increases mono- 
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Fig . 12. — z-band GCLF turnover magnitude {upper panel) and dispersion 
(lower panel) infeired from the Gaussian fits to the binned-galaxy samples 
in Fig.|6] The horizontal line in the upper panel is the same p- = -8.4 that 
characterizes the bright galaxies in the upper panel of Fig. fTTl The line in the 
lower panel is the fit of equation to the Gaussian dispersions obtained 
from fitting all 89 of our galaxies individually (cf. Fig. |9). In both panels 
the star shows values for the MW and the triangl e re presents M3 1 . Note that 
VCC 798, the bright outlier galaxy in Figs.|9]and[TT] has been excluded from 
our "binned" samples due to its excess of faint, diffuse star clusters. 

tonically from ~ 1.9 MqI'^ to Tg ~ 2.7 MqLq'. Note 
that if any of our GCs were much younger than 13 Gyr, then 
the numerical values of these mass-to-light ratios would all be 
lower (by ^ 30%-40% at an age of 8 Gyr, for example), but 
the basic constancy of T. and the systematic increase of T;, 
for redder/more metal-rich GC systems would remain. 

This has immediate implications for our plots of the GCLF 
turnover magnitudes in Figures [TT] and [12] above. In partic- 
ular, the GCs are systematically bluer, on average, in lower- 
luminosity galaxies (e.g., Peng et al. 2006a; see also Table 
|5]). Assuming that this reflects a correlation between aver- 
age cluster metallicity and galaxy luminosity (rather than one 
between cluster age and M^ gai), the typical must be some- 
what lower for GCs in faint galaxies than in bright galaxies, 
while T- is essentially the same. The fact that the Gaussian 
GC /ij scatters slightly faintward towards fainter gai should 
then reflect a modest downward scatter in the turnover mass 
scale. But in the g band, this would be balanced to at least 



GLOBULAR CLUSTER LUMINOSITY FUNCTIONS 



29 



I 





1 1 1 1 1 1 1 1 


1 1 1 1 1 
• 

/ 




/ 






/ - 






/ 






/ 






/ 






/ 






/ 






/ 






/ 






/ 




















_ — ■ 
















H h- 


— 1 \ \ \ \ 1 1 1 — 


H \ \ 1- 


- 








. 







-2 -1 






[Fe/H] 

Fig. 13. — Predicted mass-to-light ratio T (upper Panel) in solar units in 
the z (dotted line) and g (dashed line) bands, and (g— z) color {lower panel), all 
as functions of metallicity for a 1 3-Gyr old simple stellar population accord- 
ing to the PEGASE population-synthesis model (Fioc & Rocca-Volmerange 
1997) assuming a Kennicutt (1983) stellar IMF. The arrows in both panels 
indicate the minimum and maximum average GC [Fe/H] in the ACSVCS 
galaxies, as inferred from their mean (g— z) colors (see Table|5). 

some extent by the decrease in mass-to-light ratio, and [ig 
should stay more steady as a function of Mb.^s^a- 

This interpretation of the situation is confirmed in Figure 
[T4l where in the upper panel we plot the Gaussian turnover 
masses, derived from the z- and ^-band GCLF fits as just 
described, vs. parent galaxy absolute magnitude. The av- 
erage turnover magnitudes in equation ( fT9] l and the typical 
GC mass-to-light ratios in Table |5] together imply an average 
turnover mass of 



(Mto) = (2.2 ± 0.4) X 10^ 



(20) 



for the brightest ACSVCS galaxies with M^gai < -18 (here 
we have taken the absolute magnitude of the sun to be 4.5 1 in 
the z-band and 5.10 in g). The consistency in most systems 
between the turnover masses estimated from the two band- 
passes shows that, indeed, forM^^gai -18, there is an over- 
all tendency to find more GC systems with turnover masses 
somewhat below the average for giant ellipticals, by as much 
as a factor of 2 in some cases. It also implies that the de- 
pendence of GC {(g-z)) on galaxy luminosity does primarily 
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Fig. 14. — Turnover moii Mto (upper panel) and dispersion of logarithmic 
mass (lower panel) implied by our Gaussian fits to the g- and z-band GCLFs 
of individual galaxies. The turnover masses are obtained from the magnitudes 
Pg an d p- by applying the PEGASE model mass-to-light ratios summarized 
in Fig. |13| and Table|5] The dispersion in logarithmic mass is am = o-g/2.5 or 
cr-/2.5. In both plots, results from the g-band data are represented by circles, 
and results from the z-band by crosses. In the upper part of both panels we 
show the typical behaviour of error bars as a function of Mg „j. The outlying 
points at Mb ~ — 21 .2 in both panels correspond to VCC 798, a galaxy which 
has a strong excess of faint, diffuse star clusters (Peng et al. 2006b). 



reflect metallicity variations, since if GCs had very similar 
metallicities but much younger ages in fainter galaxies, the z- 
and g-band estimates of Mjo would differ by as much as the 
fitted turnover magnitudes in ^6.1.11 

For completeness, in the lower panel of Figure[T4] we show 
the g- and z-band based estimates of the Gaussian dispersion 
of logarithmic GC masses. Since cm does not depend on the 
cluster mass-to-light ratio, but is just the magnitude dispersion 
divided by 2.5, this plot is completely equivalent to Figure|9] 
Thus we have also drawn in equation (fTTI i above, multiplied 
by 0.4. 

An interesting corollary to all of this is that the reliabil- 
ity of the GCLF as a distance indicator would appear to be 
somewhat bandpass-dependent, at least when applied to sub- 



galaxies with Mb, 



aal 
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19. We have just argued that the 



near-universality of the turnover magnitude in the g-band — 
and thus in the very closely related V band — is at some level 



30 



JORDAN ET AL. 



in 



in 
O 



-+ + 



g: + 
z: o 



-22 



-20 



-If 



-16 



m 








1 












g: + 






- 


o 


+ 


z; o 


- 






+* 










- to 

+ 




+ 






in 




o 
































- 

° + 




n 

+ 
















- 




Q 













+ 










cv 


-o + 


o 

+ 


+ 


o + 


- 




- + 






+ 












+ 












o ^ 
















in 




















1 





-22 



-20 



-If 



-16 



M 



B, gal 



M 



B, gal 



Fig. 15. — Left: GCLF turnover mass Mjo, inferred from the maximum-likeli hoo d estimate of Gaussian fi, as a function of Mg.g-A for binned-galaxy samples. 
MjQ has been infen'ed from fj, using the PEGASE model, as summarized in Fig.[T3]and Table|5] Crosses indicate values of Mjq obtained from fig, while circles 
ai'e the values of A/to obtained from fj,-. (Right) GCLF turnover mass Mto versus galaxy absolute magnitude, as inferred from the fits of evolved Schechter 
function to the binned-galaxy samples (using eq. (9) and PEGASE model mass-to-light ratios). The trend of decreasing Mjq with decreasing galaxy mass is as 
in the left panel, show ing again that the choice of functional form does not affect our results. Note that VCC 798, responsible for the outlying points at bright 
Mg gal ~ -21.2 in Fig. 1 141 has been excluded from our "binned" samples due to its excess of faint, diffuse star clusters. 



the fortuitive consequence of quantitatively similar decreases 
in both the turnover mass and the typical GC mass-to-light 
ratio in smaller galaxies. At longer wavelengths, however, 
mass-to-light ratios are not so sensitive to GC metallicity, and 
variations in the turnover mass carry over more directly into 
variations in turnover magnitude. We will explore this is- 
sue in more detail in future work. However, any such prag- 
matic concerns about the precision of the GCLF peak magni- 
tude as a standard candle should not detract from the main 
point of physical interest here: although the differences in 
GCLF turnover mass that we find are real, they are never- 
theless relatively modest. While the galaxies in the ACSVCS 
range over a factor of ~ 400 in luminosity, Mto never falls 
more than ~ 30%^0% away from the (Gaussian) average of 
2.2 X 10^ Mq for the giant ellipticals. 

In the left panel of Figure[T5]we show the turnover masses 
derived from the Gaussian GCLF means for our binned- 
galaxy GC samples. This again highlights the tendency to 
slightly less massive GCLF peaks, on average, in lower- 
luminosity galaxies. In the right panel of this figure we also 
show Mto as derived from our fits of an evolved Schechter 
function to the same GCLFs (see eqs. ||9l and ifTOl ). The close 
similarity of the two graphs in Figure [15] is entirely in keep- 
ing with the slight average offset between the Gaussian and 
extended-Schechter turnover magnitudes in Figure [8] above. 
It also illustrates that our main results are not overly depen- 
dent on the particular choice of model to fit the GCLFs. 

Last, in Figure [16] we show the GC mass scales (the 
high-mass exponential cut-off) and A (interpreted as the av- 
erage mass lost per GC by evaporation) for our fits of evolved 
Schechter functions to the binned-galaxy GCLFs, as inferred 
from their magnitude equivalents and 6 in Table [3] The 
upper panel of the figure first plots M^. vs. M^ gai, using solid 
points to represent fits to GC samples selected on the basis 
of our probabilities pcc and open symbols for fits to samples 
defined only by cuts on magnitude and GC effective radius 
(see ^4.3l l. There is a clear, systematic decrease of M^ with 



decreasing galaxy luminosity. In terms of the structure of the 
mass function (eq. Q), this corresponds to a steeper fall-off 
in the frequency of GCs more massive than the turnover point. 
It is therefore equivalent to our findings in Figures [9] and [TO] 
that the Gaussian a is narrower, and the high-mass power law 
/3 steeper, for the GCLFs in fainter galaxies. As we discuss in 
^ features such as this likely reflect the initial condition of 
the GC mass distribution. Thus, if GC systems were indeed 
born with Schechter-like mass functions, it would seem that 
the "truncation" mass scale Mc was higher in larger galaxies 
right from the point of cluster formation. 

The graph of A vs. galaxy luminosity in the lower panel of 
Figure[T6]shows, first, that it is roughly comparable to (though 
slightly larger than) the GCLF turnover mass in general. This 
is certainly not unexpected, given the characteristics of the 
model itself (see the discussion in ^3.2t . In physical terms, 
though, if the model is taken at face value, the correspondence 
reflects the fundamental role that evaporation is assumed to 
play in defining any turnover point at all (see our discussion 
in |3] and the more detailed exposition of Fall & Zhang 2001). 
Beyond this, our fits imply that there is a tendency for A to 
increase as galaxy luminosity decreases, but this is not a par- 
ticularly regular trend. All in all, there appears to be a fairly 
narrow range of GC mass loss, A 2- 10 x 1O^M0, required 
to account for our GCLF observations over a large range of 
galaxy luminosity. 

Note that several of the faintest galaxy bins in Figure [16] 
have A/Mc w 2, to be compared with A/M^. ^0.1 for the 
brightest systems. This reflects once again the systematic nar- 
rowing of the GCLF, due to the steepening of dN/dM for high 
cluster masses, in fainter galaxies. 

In ^we will further discuss the variations of Mto, Mc, and 
A with galaxy luminosity, and how they relate to questions of 
GC formation and dynamical evolution. 

6.2. GCLF Turnovers in the Faintest Galaxies 
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Fig. 16. — Top: Cut-off mass scale Mc.z, inferred from our fits of evolved 
Schechter functions to the z-band GCLFs of binned-galaxy cluster samples, 
versus Mg ^.j^. Bottom: Average mass loss per globular cluster, A;, versus 
Mg gal, from fits of evolved Schechter functions to the z-band GCLFs of the 
binned-galaxy samples. In both panels, filled points are for fits to GC sam- 
ples defined by the criterion poc > 0.5 while open symbols are for samp les 
constructed using only cuts in magnitude and half-hght radius (see ji4.3i . It 
is clear that selection effects can be safely ignored when investigating broad 
trends in Mc and A. The stars in the two pan els show the values of Mc and 
A from fits to the Milky Way GCLF (eq. (l4l). The dotted fines in the lower 
panel sho w th e rou gh expe cted scaling of A vs. Mg.iiA indicated by equa- 
tions j25) and (26) in j|7.1| In order to show the scaling we have arbitrarily 
assumed that A = 2.5 X W^Mq at Mg gai = -21. 



In all of our galaxies there is evidence for the presence of 
a peak in the GCLF. Recently, van den Bergh (2006) claimed 
that the combined GCLF for a sample of local dwarf galaxies 
fainter than My.gai > -16 does not show a turnover, but con- 
tinues to increase to GC masses as low as w 1O'*M0. These 
galaxy luminosities translate to B-band gai > -15.2, which 
is essentially the magnitude limit of our ACSVCS sample. 

Even though we do not probe down to the galaxy luminosi- 
ties where van den Bergh (2006) claims a drastically differ- 
ent GCLF behavior, it is nonetheless worth noting that the 
turnover mass in our faintest galaxies is still fairly close to 
the "canonical" Mjo 2 x 1O^M0. There is no hint of any 
systematics that would cause Mto to fall to 1O''M0 or less 
in galaxies just 1 mag fainter than the smallest systems ob- 



served here (e.g., see Fig. [T4] |. It is thus likely relevant that 
the results of van den Bergh (2006) are based mostly on data 
from Sharina et al. (2005), who do not account for any poten- 
tial contamination in their lists of candidate GCs in the local 
dwarfs. Any GCLF derived from these data must therefore be 
regarded as quite uncertain, at the faint end especially. Spec- 
troscopic confirmation of the Sharina et al. GC candidates is 
required. 

6.3. Variations with Galactocentric Radius 

To achieve a fuller understanding of the GCLF, and in par- 
ticular the competing influences of cluster formation and dy- 
namical evolution on it, we would like to know how it might 
vary in form as a function of position in its parent galaxy. It 
has long been understood that the turnover of the Milky Way 
GCLF is essentially invariant with Galactocentric radius (e.g., 
Harris 2001), and multiple studies of the M87 GCLF have 
concluded that its overall shape is basically the same from the 
center of the galaxy out to several effective radii (McLaughlin, 
Harris, & Hanes 1994; Harris, Harris, & McLaughlin 1998; 
Kundu et al. 1999). Beyond this, however, little is known 
about the generic situation in most galaxies. 

For the most part, our data are not well-suited to address this 
question, due to the small field of view of the ACS. However, 
we are afforded serendipitously long baselines of galactocen- 
ti-ic radius in M87 and M49, by the inclusion in the ACSVCS 
of a number of low-luminosity galaxies that are projected 
close to each of these large galaxies. We refer to these galax- 
ies as "companions," even though they might not be physi- 
cally associated with their "hosts." The majority of the GCs 
observed in the fields of these smaller systems belong to the 
giants. While each companion does have some GCs of its 
own, their numbers will be reduced to negligible levels, com- 
pared to the M87 or M49 globulars, outside some sufficiently 
large radius in the low-luminosity galaxy. Thus, we take our 
original GC samples for the companions present in the survey 
and consider only those cluster candidates that are found more 
than 6 effective radii from the companion centers. Since the 
effective radii of the GC spatial distributions are generally w 2 
times larger than those of the underlying galaxy light (Peng 
et al. 2006, in preparation), this corresponds to excluding 
sources that are within about 3 GC-system scale radii from the 
companion centers. This should effectively eliminate w90% 
of each companion's native GCs, leaving us with fairly clean 
samples of extra M49 and M87 globulars, located tens of kpc 
away from the giant galaxy centers. 

We restrict our analysis to companions that have more 
than 50 GC candidates left after this selection. These 
are VCC 1199 (companion to M49, projected 4'.5 away); 
VCC 1 192 (M49, 4.'2); VCC 1297 (M87, 7f3); and VCC 1327 
(M87, T.5). Note that 1' = 4.8 kpc for an average distance of 
D= 16.5 Mpc to Virgo. 

In Figure [17] we show the luminosity functions and Gaus- 
sian fits for the resulting GC samples in the four fields neigh- 
bouring M87 and M49. In Table |6] we list the best-fit param- 
eters and the mean (g-z) colors and mass-to-light ratios as- 
sumed to convert the results to mass. The results are summa- 
rized in Figure [18] where we show the GCLF turnovers and 

Because the light profiles of the companion galaxies might have been 
affected by an interaction with their giant host (in the case they were physi- 
cally associated), we use the median effective radius of all VCS galaxies with 
magnitides within 0.5 mag of each companion galaxy, instead of their mea- 
sured one (the effective radii of all ACSVCS galaxies have been measured by 
Ferrarese et al. 2()06a). 
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dispersions as a function of galactocentric distance in M87 
and M49 separately. Evidently, none of the Gaussian GCLF 
parameters shows significant (> 3(j) variation over the 20-35 
kpc baselines probed. Fits of evolved Schechter functions to 
these GCLFs confirm that Mjp in particular does not change. 
As we discuss further in ^7.11 this lack of any significant ra- 
dial trend in Mjq with galactocentric distance is hard to recon- 
cile with a picture in which the GCLF turnover is determined 
solely by dynamical effects (primarily evaporation) acting on 
a universal power-law like initial cluster mass function evolv- 
ing in a fixed, time-independent galaxy potential. (In fact, if 
it varies at all, Mjq may even get slightly more massive with 
increasing radius in Fig.[l8] While we do not claim that any 
such trend is in fact detected here, it would be opposite to 
naive expectations.) 

7. DISCUSSION 

We have found interesting trends in three mass scales of 
physical interest in connection with GC luminosity functions. 

The GCLF turnover or peak mass takes a value of Mjq = 
(2.2 ± 0.4) X 10^ M0 in most bright galaxies, but shows some 
downward scatter in dwarfs fainter than Ms.gM Jc^ -18. In M87 
and M49, the data are consistent with a more or less constant 
Mto to projected galactocentric radii of 20-35 kpc. 

The higher-mass scale Mc in an evolved Schechter func- 
tion, which marks the onset of an exponential cut-off in the 
number of clusters per unit mass, grows steadily smaller in 
fainter galaxies. This drives a systematic narrowing of the 
dispersion in more traditional Gaussian fits to the GCLF, or 
equivalently a steepening of pure power-law fits to the mass 
function dN/dM at cluster masses M > Mjq. 

And the mass A in the evolved Schechter function, which 
controls the shape of the low-mass end of the GC mass distri- 
bution and is instrumental in setting Mjq, varies by factors of 
a few — although not entirely monotonically — as a function of 
galaxy luminosity. 

We now discuss these results in terms of their implications 
for GC formation and dynamical evolution. We begin by fo- 
cusing on A in the evolved Schechter function, which, in the 
context of Fall & Zhang's (2001) dynamical theory for the 
GCLF, is meant to measure the average amount of mass lost 
per globular cluster in a galaxy, over a Hubble time of evolu- 
tion. We then move on to Mc and Mto, asking specifically to 
what extent the observed variations in these high-mass char- 
acteristics of the GCLF might be caused by dynamical friction 
rather than initial conditions. 

7.1. Evaporation and the Low-Mass Side of the GCLF 

The defining feature of the evolved Schechter function in 
equation (jTjl — which we have found to fit the GC mass distri- 
butions of galaxies in the ACSVCS just as well as the tra- 
ditional, but ad hoc, lognormal form — is the flat shape of 
dN/dM in the limit of low masses. This asymptotic flat- 
ness always follows naturally from a time-independent rate 
of cluster mass loss, regardless of the assumed initial form 
of dN/dMo (Fan & Zhang 2001, and EI2]above). The ex- 
act values of the average cumulative mass losses per GC for 
the galaxies in our sample are, however, more specific to the 
assumption that dN/dMo oc Mq^ exp(-M()/Mc) — a form cho- 
sen to match the observed mass functions of young massive 
clusters in local mergers and starbursts. 

It is worth noting that, even though the average mass loss A 
in an evolved Schechter function is key to setting the GCLF 
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Fig . 17. — Histograms and Gaussian fits for the GCLFs of GCs in tlie field 
of view of four companions of M87 (=VCC 1316) and M49 (=VCC 1226) 
that lie farther away than 6 Re from the companion galaxy, where Re is the 
effective radius of the companion and is determined as described in the text. 
For each field we present the z-band and g-band GCLFs side by side. The 
VCC name of the companion galaxy is indicated in the upper left corner of 
the left panel, where we also indicate the total number N of sources with 
Pac > 0.5 and the bin-width h used when constructing the histograms. In 
each panel we show the best fitting model (solid black curve), the intrinsic 
Gaussian component (dashed curve), the Gaussian component multiplied by 
the completeness fraction (dotted curve), and a kernel-density estimate of the 
expected contamination in the sample (solid gray curve). The solid black 
curve is the sum of the solid gray and dotted curves. Details of the fits are 
given in Table[6] 

turnover mass, Mjq does not vary as much or as systemati- 
cally as A does in the ACSVCS sample (cf. Figures [TS] and 
[T6] ). This is because the value of the upper-mass cut-off Mc 
also influences Mjq (see ^3.21) . and Mc varies in such a way 
as to largely counteract the variation of A, keeping Mto more 
steady as a function of M^ gai. 

Since Mto is observed to be so nearly constant indepen- 
dently of any functional fitting — at least in large galaxies — 
this balance between variations in A and Mc might be viewed 
simply as a necessary condition to make evolved Schechter 
functions match the data at all. But more interesting is that 
if the physical arguments behind the fitting function are close 
to correct, our results imply that the near-universality of the 
GCLF turnover in bright galaxies (M^.gai ^ -18) is in some 



GLOBULAR CLUSTER LUMINOSITY FUNCTIONS 



33 



VCC1226 = M49 



VCC1316 = M87 



VCC1226 = M49 



VCC1316 = M87 



(kpc) 



20 30 
R , (kpc) 



Fig. 18. — (Left) GCLF turnover mass Mjo as a function of projected galactocentric distance Rgc in M49=VCC 1226 (top) and M87=VCC 1316 (bottom). 
Filled symbols refer to z-band measurements of Mjq while open symbols refer to g-band ones. (Right) Gaussian dispersion of logarithmic cluster masses, ctm, 
as a function of projected galactocentric distance Rgc in M49=VCC1226 (top) and M87=VCC1316 (bottom). Filled symbols refer to z-band measurements while 
open symbols refer to g-band ones. In the M49 panels, the leftmost pairs of points refer to the GCLF parameters derived from the central ACS pointings in the 
main body of the galaxy; the next pair out coirespond to the companion VCC 1192; and the rightmost pair correspond to the companion VCC 1199. In the 
M87 panels, the leftmost pairs of points refer to the GCLF parameters derived from the central ACS pointings in the main body of the galaxy; the next pair out 
correspond to the companion VCC 1297; and the rightmost pair con'espond to the companion VCC 1327. 



sense a coincidence resulting from steeper initial dN/dM 
(with lower M^) in fainter systems being eroded by faster 
mass-loss rates (yielding larger A). 

As we discussed in ^3.21 some amount of cluster mass loss 
may result from tidal shocks, but we expect that in general 
the largest part comes from two-body relaxation and evapora- 
tion, at a rate determined by the mean cluster density inside 

its half-mass radius: /iev oc p'/^. This basic dependence holds 
independently of any host galaxy properties, so if the cluster 
evaporation rate varies systematically as a function of Mb.^a, 
it presumably reflects systematics in the typical pi, of the clus- 
ter systems. Then, if GCs are tidally limited, such that their 
average densities are determined by the galaxy density inside 
their orbits (e.g.. King 1962), variations in their character- 
istic ph should correspond in some way to variations in the 
host-galaxy densities. The easiest way to quantify any such 
connection is to assume a spherical, time-independent galaxy 
potential with a simple analytical form. Thus, in their mod- 
els of the Milky Way GC system. Fall & Zhang (2001) relate 
the ph of individual clusters to their orbital pericenters rp in 
a logarithmic potential with a fixed circular speed, Ve, so that 
Ph OC Pgaiirp) (X jr-p. We address the validity of these partic- 
ular (strong) assumptions about the host galaxy below; but for 
the moment we follow Fall & Zhang and most other authors 
(e.g., Vesperini 2000, 2001; Baumgardt & Makino 2003) in 
making them. What do our fitted A values for the ACSVCS 
galaxies then imply for the distribution of GC densities and 
pericenters in these systems? 

The evaporation rate of a cluster with observable, projected 
half-mass radius Rk depends on the density ph -~ 
roughly as 



galaxy potential.^' However, if p/, is taken to be set by a well 
defined Vp in a steady-state, singular isothermal sphere, then 
we also have (from eqs. [4] and [15] of Fall & Zhang 2001) 



Pev(theo)~2.9x lO'^MoGyr"' (r^/kpc)" 
X [l-ln(rp/r,)]'/' . 



(K/220kms"') 
(22) 



In the last term on the right-hand side, which is derived by In- 
nanen, Harris, & Webbink (1983), is the radius of a circular 
orbit with the same energy as an arbitrary orbit with Vp < re- 

Now, for the Milky Way, recall from ^33] (eq. |[T4l ), that 
we estimate 



A(MW) = (2.5 ± 0.5) X 10^ Mq 



(23) 



from our fit of an evolved Schechter function to the GCLF. For 
a GC age of 13 Gyr, this implies a mass-loss rate (averaged 
over the distribution of cluster pi, or, given the assumptions 
behind eq. 1221 . over all cluster orbits) of 



(Aiev)(fit) = 



A(MW) 
13 Gyr 



:(1.9±0.4)x lO^MeGyr"' 



(24) 



3M/i8TrRl) 



Pev(theo) ~ 345 Mq Gyr 



Ph 



Mopc- 



1/2 



(21) 



which again is independent of any assumptions on the host- 



Comparing equation ( l24b to equation ( 1211 1 implies an average 
(ph) ~ (3000 ± 600) Mq pc-3 for GCs in the Milky Way. This 
average falls towards the upper end of the range of cluster p/, 
observed today, but it is within a factor of ~ 2-3 of the mean 
(e.g., see the data in Harris 1996). Equation (l22l) further sug- 
gests an average pericenter of (rp) w 2 kpc. This is roughly 
the same answer found by Fall & Zhang (2001; see their Fig- 
ure 13), which shows that an evolved Schechter function is 

^' Equation )211 follows from /Xcv = 0.045M/fih (Henon 1961; see also 
Fall & Zhang 2001), where the half-mass relaxation time ^.i, is given by equa- 
tion (8-72) of Binney & Tremaine (1987) with (1) an average stellar mass of 
0.1 Mq and a Coulomb logarithm In A = 12 assumed constant in time, as in 
Fall & Zhang, and (2) a generic proportionality, Rj, ~ 0.75 r;,, between the 
projected half-mass radius R/, and its unprojected counterpart r/, (e.g., Spitzer 
1987). 
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a reasonable analytical approximation to their full numerical 
theory. While such an {rp) is slightly small — just as (/?/,) is 
slightly high — compared to more direct pericenter estimates 
for Galactic globulars (cf. Innanen et al. 1983; van den Bergh 
1995), it is again within the range of standard values. 

It is not at all obvious a priori that average cluster densi- 
ties and pericenters inferred strictly from fits to the Galac- 
tic GC mass function should agree to within factors of 2 
or 3 with values estimated by independent methods. The 
fact that they do is an encouraging sign for the basic picture 
of evaporation-dominated GCLF evolution. Some residual 
corrections — downward in "predicted" (p/,) and up in (r^) — 
are evidently required, but at a level that plausibly could come 
from straightforward refinements in the various steps lead- 
ing to equations ( 1211 1 and ( l22b . For example, there is some 
room for adjustment of the exact theoretical coefficients for 

the evaporation rate /^ev « M/f rh c>c p,/ and the pericenter 

rp (X (see, e.g., the discussions in Fall & Zhang 2001). 
In addition, we have neglected here any additional mass loss 
caused by tidal shocks, and we have adopted the idealization 
of a spherical and time-invariant Galactic potential. 

To bring the ACSVCS data into this discussion, we focus 
on the basic pattern of variation in A as a function of galaxy 
luminosity, shown in the lower panel of Figure[T6] First, A in- 
creases slightly from the brightest Mfi gai ~ -21 .5 to the fainter 
Mb.hM — -18. The uncertainties and scatter in A are large, 
but the mean increase is perhaps a factor of «2-5. Then, at 
fainter Mfi.gai > -18, A holds more constant or even decreases 
again, possibly by as much as a factor of w2-3 by the limiting 
Ms.g-di — -16 of the survey. 

If evaporation is responsible for these variations, then we 
should expect them to be mirrored in the behavior of the aver- 
age GC half-mass radius as a function of galaxy luminosity: 
from equation (l2Tl i. (/iev) oc (p/,)''^^ oc {Rh)~^/^, and by defi- 
nition A cx (/Xev) for coeval clusters. Globulars in Virgo are 
marginally resolved with the ACS, and Jordan et al. (2005) 
have fit PSF-convolved King (1966) models to estimate in- 
trinsic Ri, values for individual sources (selected as described 
in §2 of Jordan et al. 2005) in most of the galaxies that we 
have dealt with here. The behavior of mean {Rk) versus gai 
is shown in Figure 5 of Jordan et al. 

A detailed comparison of (7?/,) and A is not straightfor- 
ward, since these quantities were estimated separately for 
GC samples defined differently by Jordan et al. than in this 
paper Nevertheless, it is interesting that (/?/,) can be de- 
scribed as decreasing towards brighter galaxy luminosity in 
the range -21.5 < M^aai J; -18, where A increases, and then 
turning around to increase somewhat at fainter M^ gai > -18, 
where A decreases again. The changes in {Ri,) are — as we 
would expect — smaller and less clear than those in A, but it 
is just plausible that there is a net decrease of ~ 35% from 
■^B.gai = -21.5 to Mb gal = -18 and a slightly larger increase 
from Mfiaai = -18 to M^^gai = -16. This would be consistent 
with the shallowest trends able to fit A versus Mb gai in Figure 

m 

We cannot use equation ( |22] | to relate A to typical GC peri- 
centers and average galaxy densities on a case-by-case basis 
in the ACSVCS sample as in the Milky Way, since Vc obser- 
vations are not available for all systems. However, scaling 
relations can be used to some effect here. Large early-type 
galaxies with M^ gai ;$ -18 generally obey Vc = V2a oc L^^-^^ 
(e.g., Faber & Jackson 1976), (M/L)gai cx L°;,f at optical 



wavelengths (van der Marel 1991; Cappellari et al. 2006), 
and thus R^s (x igaJ"^ by the virial theorem (see also Ha§e- 
gan et al. 2005). Average mass densities therefore increase 
towards lower Lgai, such that equation (l22l i implies 

A (X {R.ff/{rp)) (X (/?eff/(r,» LgaV^"'"' (25) 

for bright galaxies. The situation is somewhat different for 
fainter M^gai ^ -18. For Coma Cluster galaxies in this 

regime, Matkovic & Guzman (2005) find Vc = \/2cr cx Lgai^° \ 
while the data in Graham & Guzman (2003) suggest Rf^s cx 
^gai"^- these systems are representative of those in Virgo, 
then their average densities decrease towards lower Lgai, and 
equation ( l22b leads to 

A cx {R,s/{rp)) L^^f (26) 

for faint dwarfs. 

The major unknown in equations ( |25] ) and ( |26] l is the ratio 
of galaxy 7?eff to GC (r^), and how it might or might not vary 
systematically as a function of galaxy luminosity. If the ratio 
is constant for all systems, then the dotted lines drawn in the 
lower panel of Figure [16] show the expected variation of the 
mass loss A versus gai. These lines are normalized to make 
A = 2.5 X 10^ M0 at M^.gai = -21 and to make the bright- and 
faint-galaxy scalings meet at M^ gai = -18. The net increase of 
^gai^ fromMfi gai = -21.5 toM^ gai = -18 is a factor of about 5, 
while the decrease of Lg^^-]^ from M^ gai = -18 to M^.gai = -16 
is a factor of approximately 2. 

These changes may be somewhat greater than suggested by 
the actual fitted estimates of A. Moreover, an increase of A by 
a factor of 5 between M^ gai = -21.5 and M^ gai = -18 would 
imply a decrease in (/?/,) by a factor of 5^^^ « 3, which is 
larger than the measurements of Jordan et al. (2005) support. 
However, this is clearly not an order-of-magnitude problem. 
It could easily be alleviated if the galaxy total mass distribu- 
tions are not isothermal spheres, or if Rss/{rp) depends even 
weakly on galaxy luminosity, or if uncertainties and scatter in 
the galaxy scalings result in small deviations from the nomi- 
nal exponents on Lgai in equations ( l25T l and ( l26b . Tidal shocks 
may also contribute differently to the net A in different galax- 
ies, a complication that we have entirely ignored. Again, then, 
it is encouraging that these crude relations come as close as 
they do to explaining the systematics in a cluster mass-loss 
parameter inferred only from the GCLF — accounting in par- 
ticular for the change in dependence of A on galaxy luminos- 
ity around Mfi.gai ~ -18. 

Obviously, more rigorous and detailed analyses of individ- 
ual galaxies are required to really make (or break) the case in 
general that the overall form of an evolved Schechter function 
for the GC mass function, and the parameter A especially, can 
be interpreted physically and self-consistently as the result of 
evolution from an initial GC dN/dMo cx My^ with individual 
cluster mass-loss rates that are constant in time. From our dis- 
cussion here, it does seem that this "literal" view of the simple 
fits to the Milky Way and ACSVCS GCLFs is at least broadly 
compatible with observations of the cluster densities or radii 
in these galaxies and with the trends in A vs. Lgai, if evapora- 
tion is the main disruptive process for clusters as massive as 
Mjo^lx IO^Mq. 

Difficulties do arise, however, when considering the addi- 
tional constraint that the GCLF is invariant over wide ranges 
of galactocentric radius and GC density in the Milky Way and 
other large galaxies. As described above, application of equa- 
tion ( |22] ) to the global Galactic GCLF ultimately implies an 
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average GC pericenter of (r,,) ~ 2 kpc, corresponding to about 
half the effective radius of the bulge. Similarly, our normal- 
ization of equation (l25T l in Figure [16] implies (r^) < 0.5/?eff 
for the brightest early-type galaxies in Virgo. But observa- 
tionally, the GCLF turnover Mjo (and thus A) has the same, 
global value for clusters currently found out to at least 10- 
15 effective radii in the Milky Way (e.g., Harris 2001) and at 
least ~ 4Reff in M87 and M49 (E3- This can only be consis- 
tent with evaporation-dominated depletion of an intially steep 
GC dN /dM() cx My- at low masses, and with the additional 
assumption that the mass loss A (x r^', if cluster orbits are 
systematically much more elongated at larger galactocentric 
radius in all these systems. 

In fact, for the Milky Way and M87 respectively. Fall & 
Zhang (2001) and Vesperini et al. (2003) have shown that fol- 
lowing this chain of logic leads to the conclusion that globu- 
lars should initially have been on predominantly radial orbits 
outside about one effective radius in each galaxy. On the other 
hand, the present GC velocity distributions in the Galaxy, in 
M87, and in M49 are all essentially isotropic — implying or- 
bits with typical axis ratios of only ro/r^ ~ 3 — out to the same 
spatial scales of several /?eff, over which the observed GCLF 
is invariant (see, e.g., Dinescu, Girard, & van Altena 1999; 
Cote et al. 2001, 2003). Fall & Zhang suggest that this dif- 
ference between (presumed) initial and (observed) present or- 
bital properties might be explained by preferential depletion 
of GCs on the most radial orbits. But while the idea remains 
to be tested in detail for the Milky Way, Vesperini et al. (2003) 
show that — again if the galaxy potential is spherical and time- 
independent — it does not suffice to account quantitatively for 
the combined GCLF and kinematics data in M87. 

Related to this is the average density, (p;,) ~ 3000 MqPC"^, 
implied by the more general equation (l2ll i and the required to- 
tal A for Galactic globulars. A similar (p/,) is also suggested 
for GCs in the brightest Virgo galaxies by the A values in Fig- 
ure [16] As we mentioned above, such densities are observed 
for real clusters; but there is a broad distribution of ph, with 
an average slightly lower than 3000 Mq pc~^ and a long tail 
to much smaller values of < 100 Mq^c~^. More generally, 
the GCs in most large galaxies have half-mass radii that are 
largely uncorrected with cluster mass (e.g., van den Bergh, 
Morbey, & Pazder 1991; Jordan et al. 2005, and references 
therein), so that pi, apparently always ranges over more than 
two orders of magnitude. When p/, < 100 Mopc"^, the total 
evaporative mass loss per cluster over 13 Gyr is <5x 1O'*M0, 
well below the typical average A and global Mjo for entire 
GC systems. In the Milky Way at least, the large majority of 
such low-density GCs are found at Galactocentric distances 
Tgc > 10 kpc, so in a sense the problem is bound up with the 
weak radial variation of the GCLF. 

These points are important, and they need to be resolved, 
but they should not be taken as disproof of the idea that 
long-term dynamical evolution alone might explain the dif- 
ference between the mass functions of old GCs and young 
massive clusters. Ultimately, the near-flatness of dN jdM at 
low masses, which is clearly seen in the Milky Way and is en- 
tirely consistent with all of our Virgo GCLFs, only demands 
that cluster masses decrease linearly in time if the dynamical- 
evolution hypothesis is correct at all (see Fall & Zhang 2001, 
and 33.21 above). It is not absolutely necessary that evap- 
oration account for the full mass-loss rate of every cluster, 
even though our discussion here has focused on exploring this 
possibility (and shown that it does come remarkably close. 



to within factors of 2-3 for the most part). For example, 
globulars in the extreme low-density tails of ph distributions, 
mentioned just above, might be much more strongly — and 
differently — affected by tidal shocks than any previous GCLF 
calculations have allowed. Such shock-dominated evolution 
could still lead to a cons tant m ass-loss rate of its own (see 
Dehnen et al. 2004, and 33.2.21 above), which would add di- 
rectly to /iev without otherwise changing any of the main ar- 
guments here. 

In more specific terms, the radial invariance of the GCLF 
might ultimately be explained by modifying a single ancillary 
assumption in the current dynamical-evolution models rather 
than discarding the idea altogether. It is the notion of spherical 
and steady-state galaxy potentials that prompts Fall & Zhang 
(2001), Vesperini et al. (2003), and almost all other authors to 
use equations (1211 1 and (l22l l to tie cluster densities to orbital 
pericenters in these analyses. But, as Fall & Zhang them- 
selves point out, this is of course an extreme simplification 
for galaxies that grow through hierarchical merging. 

Fall & Zhang suggest, for example, that a major merger 
could obviate the need for extremely radial orbits to distribute 
clusters with high mean densities, fixed at small and well de- 
fined pericenters, over large volumes in a galaxy. Instead, a 
merger may efficiently mix two globular cluster systems spa- 
tially and isotropize their velocity distribution. This could 
then work to weaken any radial gradients in the mass loss 
A and the GCLF turnover mass, which might have resulted 
from realistic orbital distributions and Pev{rp) dependence like 
equation (l22l i in the progenitor galaxies. 

In addition to this, multiple minor mergers — which are per- 
haps more relevant than major mergers for a galaxy like our 
own — should steadily bring in globulars formed with densi- 
ties and evaporation rates unrelated, at least initially, to their 
new orbits in the main galaxy, making the use of equation (i22l i 
less than straightforward. In fact, any use of it at all could be 
questionable in this case, since all clusters would constantly 
be sampling new pericenters in an evolving potential. Again, 
then, weak spatial variations in A and Mjo need not imply 
highly radial GC orbits. Prieto & Gnedin (2006) have recently 
simulated the evolution of the GCLF during the hierarchical 
growth of a Milky Way-sized galaxy. Starting from an ini- 
tial cluster mass function dN /dMo cx M^^, which is re-shaped 
primarily by evaporation — ^but abandoning equation (l22l i and 
instead adopting evaporation rates from GC densities fixed in- 
dependently of their orbits — they find that it is possible (even 
without a recent major merger) to produce a final GC system 
with an isotropic velocity distribution and a radially invariant 
GCLF similar to the observed Galactic distribution. 

A caveat is that the hierarchical-growth simulations most 
favored by Prieto & Gnedin (2006) are ones in which they as- 
sume that all globular clusters have a common mean density 
inside Ri, (just one that is not set by any orbital pericenter). 
This is still incompatible with the wide range of ph observed 
for the GCs in many galaxies, and it is furthermore not ob- 
vious how the cumulative mass loss A cx {phY^^ should then 
vary as a function of galaxy luminosity. On the other hand, 
Prieto & Gnedin have also run some models allowing for an 
initial spread of GC densities followed by evaporation at con- 
stant pij. This is at least more reminiscent of real p/, distribu- 
tions, and it still produces a GCLF that is not too drastically 
different from the Galactic one. Clearly, more work is re- 
quired to clarify the dynamical evolution of initial power-law 
GC mass functions in time-dependent galaxy potentials, with 
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the totality of relevant observational constraints taken into ac- 
count: a flat dN/dM at low masses; a weak or absent correla- 
tion between GC radii and masses; radially invariant GCLFs; 
currently isotropic velocity distributions; and mass losses A 
that vary with galaxy luminosity as in Figure [T6l 

Should all efforts along these lines fail to explain the com- 
bined data, the only option left would seem to be that a peak 
in the GCLF was established much earlier, by processes more 
related to cluster formation. One possible scenario has been 
proposed by Vesperini & Zepf (2003). They suggest that low- 
mass globulars were initially less concentrated (with a larger 
ratio of half-mass to tidal radius) than high-mass clusters. The 
inevitable expansion of all clusters following mass loss driven 
by stellar evolution would then cause many low-mass clus- 
ters preferentially to overflow their tidal radii, leading ulti- 
mately to fast disruption times of a few hundred Myr or less 
(Chernoff & Weinberg 1990). This may turn an initial power- 
law dN/dMo at low masses into a roughly flat-topped or even 
lognormal distribution, with Mjo near its current value, very 
early on. Weaker long-term evaporation (i.e., lower cluster 
densities or larger and more variable pericenters) could then 
suffice to explain the residual difference between the initial, 
steep mass function and the final, observed one, even in a 
static galaxy potential. 

Observations of the young massive clusters in the Antennae 
galaxies already imply that early disruption is independent of 
cluster mass, at least for clusters more massive than several 
IQ'^Mq and younger than ~ 10*^ yr (Zhang & Fall 1999; Fall 
et al. 2005). Thus, if the disruption mechanism of Vesperini 
& Zepf (2003) is to work, the mass-selective aspect of it ap- 
parently must be restricted to timescales of lO'^-lO'' yr or so. 
In any case, the success of this or any similar picture further 
relies on an appropriately tuned mass dependence in some key 
GC property being built into cluster systems essentially as an 
initial condition; but this still requires explanation in itself. 




log (M/Mq) 

Fig. 19. — A simple model for the effects of dynamical friction on the GC 
mass function, illustrated in terms of the distribution ilN /d log M, which is 
directly proportional to the GCLF. The solid curve shows this version of an 
evolved Schechter function with Mc = 3.0 X 10^ Mq and A = 2.6 X 10^ Mq 
(hence Mjq = 2.2 X 10^ Mq), appropriate for a giant elliptical. This is as- 
sumed to be the mass function on which dynamical friction operates. The 
mon oton ically decreasing dotted curve shows the function S(M.t = 13 Gyr) 
(eq. 1301 ). calculated for a galaxy with Mg gai = —21.75 as described in the 
text. The dotted bell-shaped curve is the product of this S times the solid 
curve; it illustrates the cumulative effect of dynamical friction on the GCLF 
in very massive galaxies. The monotonically decreasing dashed curve shows 
the function S(M,t =13 Gyr) as calculated for a galaxy with Mg gai = —15.75. 
The dashed bell-shaped curve is the product of this S times the solid curve, 
illustrating the net effect of dynamical friction on the GCLF in very faint 
galaxies. The arrows indicate the position of the final turnover mass for each 
of the resulting mass functions. There is a sHght decrease in Mxo for the low- 
luminosity galaxy as a consequence of dynam ical frictio n, bu t not enough to 
account fully for the observed behavior in Fig. II 5 l or Fig. 1 121 



7.2. Dynamical Friction and the High-Mass Side of the 

GCLF 

At GC mass scales M > A, dynamical friction can in some 
cases become more important than evaporation or shocks as a 
cluster destruction mechanism. A point mass M originally on 
a circular orbit of radius r in a galaxy with a total-mass distri- 
bution following a singular isothermal sphere will spiral in to 
the galaxy center within a time (Binney & Tremaine 1987) 



Tdf : 



5.9 Gyr / r \ 



(lnA)/10 Vkpcy V220kms 



Vr 



10* 
M 



, (27) 



where Vc is the galaxy's circular speed and InA ~ 10 is the 
usual Coulomb logarithm. 

It is clear from equation (|27] i that dynamical friction can- 
not be a major factor in deciding the evolution of all but the 
very most massive tip of the GCLF in and brighter galax- 
ies with Vc > 200 km s"'. However, the scaling Tdf cx Vc im- 
plies that the relevance of dynamical friction can increase sig- 
nificantly for lower luminosity galaxies (e.g., Hernandez & 
Gilmore 1998; Lotz et al. 2001). It is then reasonable to ask 
whether a stronger depletion of massive GCs in dwarf galax- 
ies might be able to explain the systematic decrease of Mc ver- 
sus Mb, gal in our fits of evolved Schechter functions for these 
systems, and possibly even the slight decrease in average Mjo 
towards the faintest Ms^gai- 

We do not attempt here to find a definitive answer to this 
question, but only an indication of the ability of dynamical 



friction to produce the observed trends. One particular sub- 
tlety is that the expression for Tdf in equation dZTl i does not al- 
low for clusters to evaporate. But a steadily decreasing cluster 
mass will lead to a longer total dynamical-friction timescale. 
We deal with this complication in the simplest way possible: 
the timescale Tdf for a cluster with initial mass Mq and present 
mass M = (Mo - A) is approximated by evaluating equation 
( |27] | at the average mass, (M+ A/2). 

Let us denote by ^'(M, f ) the GC mass function that would 
be obtained after a time f of GC evolution in the absence of 
any dynamical friction. The effects of dynamical friction are 
easily accounted for by subtracting from "if all clusters with 
instantaneous masses M such that 

(M+A/2)>M^i„(r,f), (28) 
where Mmin follows from equation (l27b by setting Tdf < f : 

4.5x10^Mq ^ 13Gyr ^ ^ MV 
M^Ur,t)- ^j^^^^^Q ^ J^^^^J l,220kms-' 

(29) 

The net, "global" GC mass function (averaged over all GC 
orbits, or galactocentric radii) at any time t is thus dN/dM = 
S(M, t) X $(M, t) where 



5(M,f) = 



/o°° PGc(r)//[M„in(r, t) - (M+ A/2)] A^r^dr 
Pcc(r)4TTr^dr 



(30) 
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Here pccir) is the space density of GCs (assumed to be inde- 
pendent of cluster mass) and H is the Heaviside step function: 
H{x) = 1 for jc > and H{x) = for jc < 0. 

This raises further points to be dealt with in more care- 
ful calculations along these lines. First, dynamical friction 
will clearly affect also the spatial distribution of GCs, so that 
PGcif) will have a dependence on time, which we ignore. Sec- 
ond, the effects of dynamical friction could introduce some 
dependence on galactocentric position into the GC mass func- 
tion, which in a complete treatment would be contrasted with 
observational limits on any such variations. Third, changing 
the assumed galaxy potential could significantly affect the de- 
rived Tdf (e.g., Hernandez & Gilmore 1998; Read et al. 2006), 
as could relaxing the unrealistic assumption of strictly circu- 
lar orbits (e.g., Pesce, Capuzzo-Dolcetta, & Vietri 1992; van 
den Bosch et al. 1999). Finally, we do not take into account 
the fact that the ACS has a fixed field of view, and thus we 
are not always observing truly globally averaged GCLFs — 
although this point is relevant mainly for the most massive 
galaxies, where the effects of dynamical friction are expected 
to be negligible in any case. 

These issues notwithstanding, we proceed to estimate the 
effects of dynamical friction by evaluating S(M,t) as writ- 
ten in equation ( l30b . We assume that the "friction-free" 
^'(M, f ) at the present day is well described by the GCLF of 
bright ellipticals, where dynamical friction is negligible, and 
is therefore given by equation (|7]l with A = 2.6 x 1O^M0 and 
Mc = 3 X 1O*'M0 (see Figure [16] and Tabled. To obtain the 
final dN/dM including dynamical friction, we then multiply 
this by the function S(M,t = 13Gyr). In doing so, we always 
take the slowly varying Coulomb logarithm in equation ( |29] l 
to be lnA= 10. 

We assume that for giant galaxies with Mg < -18 we have 
Vc oca (X Lg^'j^ (Faber & Jackson 1976), with a zeropoint cho- 
sen to give Vc = 484 km s"' at Mb = -21 .75, based on the ve- 
locity dispersion of M87 (Bender, Saglia & Gerhard 1994). 
We impose a change in this scaling at Mb > -18, so that 
dwarfs follow K cx cr cx L°J^^"■^ (Matkovic & Guzman 2005; 
cf. ^7.1l above). We can then find M^ninif, t) from equation ( l29l l 
for any GC in any galaxy. 

To specify the spatial distribution of GCs and calculate 
5(M, f = 13 Gyr), we estimate the galaxy's effective radius /?eff 
using the data from Ferrarese et al. (2006a); then we assume 
that the effective radius of the GC system is just twice R^a 
(Peng et al. 2006, in preparation). Finally, we assume that 
PGcif) is given by the density profile of Prugniel & Simien 
(1997; see also Terzic & Graham 2005), which is an analyti- 
cal approximation to the deprojection of a Sersic profile (/?'/" 
law), and we let the Sersic index n be determined by M^ gai as 
per equation (25) of Ferrarese et al. (2006a). 

The results of the calculations for two representative galaxy 
magnitudes, M^gai = -21.75 and Mb^a = -15.75, are il- 
lustrated in Figure [19] The figure shows both S(M,t = 
13 Gyr) (the monotonically decreasing curves) and the func- 
tion dN/dlogM (proportional to the GCLF and given by the 
peaked curves) that follows from dynamical friction acting 
on the assumed evolved Schechter function. The resulting 
turnover mass scales are indicated with arrows, which show 
that the stronger dynamical friction in the fainter galaxy leads 
to a slightly lower turnover mass scale. 

We show the behavior of Mjo as a function of Mg.gai in gen- 
eral, in the upper panel of Figure [20] (circles connected by a 
solid Une) and contrast it with the observed variation in our 
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Fig. 20. — (Upper panel) GCLF turnover mass Mto.c inferred from 
evolved Schechter fits to z-band data for the binned-galaxy GC samples, ver- 
sus Mg (data are the same as the circles in Fig. 115) . The open symbols 
connected by a solid line are the predicted change in Mjo due to the increas- 
ing efficiency of dynamical friction in fainter galaxies; see text. Dynamical 
friction might account for part of the observed trend, but probably not all of 
it. (Lower panel) Evolved Schechter function cut-off mass Mc.z versus Mg ^.^i 
for binned-galaxy z-band samples. The open symbols connected by a solid 
line are again our predictions for the change in Mc due to stronger dynami- 
cal friction in smaller galaxies. Dynamical friction is not able to explain the 
observed behavior of Mc as a function of galaxy luminosity. 

binned-galaxy GC samples (Figure [TSjl- The predicted Mjo 
varies quite slowly with Mb. gai, but it ultimately decreases by 
^ 10% from our assumed 2.2 x 10^ M0 in the brightest galax- 
ies. This is comparable to the observed decrease of 30% 
in Mto- Thus, dynamical friction may be responsible for 
some part of the the slow change in GCLF turnover mass with 
galaxy magnitude. 

In the lower panel of Figure [20] we show (again, as open 
circles connected by a solid line) the Mc values inferred by 
fitting evolved Schechter functions to our model GC mass 
functions after calculating the effects of dynamical friction. 
Evidently, we can expect dynamical friction to cause perhaps 
a ^ 30%-40% decrease in the value of from the bright- 
est to the faintest galaxies; but this is altogether too little to 
account for the factor of ~6-7 decrease we actually observe. 
Similarly, if we fit power laws to our dynamical-friction mass 
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functions in the range 3 x 10^ Mq < M < 2 X 10^ Mq, we 
obtain rather constant powers /3 ~ 1.7-1.8 for galaxy mag- 
nitudes -21.75 < Mb gal < -15.75, which is far from being 
able to explain the observational situation in Figure [TO] above. 
We conclude that dynamical friction cannot account for more 
than a small fraction of the observed steepening of the globu- 
lar cluster mass function above the GCLF turnover. 

These results are essentially in agreement with those of Ves- 
perini (2000), who models the effects of evaporation and dy- 
namical friction on the GCLF, and predicts only slight de- 
creases in the mean (log M) and the Gaussian dispersion am 
as galaxy luminosity decreases (see his Figure 6), at levels 
much smaller than those seen in our data (e.g.. Figure O. 
Thus, although we have emphasized the highly simplified na- 
ture of our calculations, it nevertheless appears that galaxy-to- 
galaxy systematics in the cluster formation processes, rather 
than dynamical evolution, must be largely responsible for the 
observed variation in the detailed form of the GCLF at high 
masses. 

7.2.1. Initial Conditions 

It seems inevitable from the discussion above that the 
observed steepening or narrowing of the GCLF above the 
turnover point Mto ^ 2 x IQ^Mq in fainter galaxies — 
whether this is expressed in terms of smaller Schechter- 
function mass scales Mc or steeper power-law indices (3 
or narrower Gaussian dispersions ctm — must reflect non- 
universal initial conditions in the cluster mass distribution, 
and therefore some fundamental aspect of the star-formation 
process. 

Observationally, it is known that the luminosity of the 
brightest young star cluster in a star-forming galaxy scales 
with the global star formation rate (Billet, Hunter & 
Elmegreen 2002; Larsen 2002). There has been some discus- 
sion as to whether this is just a size-of-sample effect (if more 
clusters are formed, it is statistically more likely to achieve 
higher masses by random sampling of an underlying mass 
distribution that might still be universal) or indicative of a 
real, physical limit to the initial cluster mass function (Larsen 
2002; Weidner, Ki-oupa & Larsen 2004). 

Gieles et al. (2006a, b) argue that there is a physical upper 
limit, Mnjax, to cluster masses in each of NGC 6946, M5 1 , and 
the Antennae galaxies (though see Whitmore, Chandar & Fall 
2006 for a differing view). The number of clusters found with 
M > Mn,ax falls rapidly to zero in all three cases, but the value 
of the upper limit is found by Gieles et al. to vary between 
the galaxies, from M^ax — 4-10 x IO^Mq. QuaUtatively, a 
parameter like M^ax can be identified with Mc in a Schechter- 
function description of (initial) GC mass functions. Quanti- 
tatively, the range of M^ax claimed by Gieles et al. for their 
young systems is very similar indeed to our fitted values 
for the old GCs in early-type Virgo galaxies (see Figure [T6ll. 

It will be interesting to explore this possible connection be- 
tween globulars and young massive clusters in more detail. 
Possibly one route to take is suggested by the theory of the 
GCLF developed by Harris & Pudritz (1994), in which a dis- 
tribution of cluster masses is built up by collisions between 
gaseous protoclusters. McLaughlin & Pudritz (1996) suggest 
that the total time required to produce very high-mass clus- 
ters may be longer for galaxies in lower-density environments, 
and this could perhaps be related to our finding of a cut-off 
at lower (in our current notation) for the initial GC mass 
functions at fainter M^aai. If these types of ideas can be gen- 
eraUzed, then both our GCLF observations and the possible 



existence of an upper mass "limit" in young cluster systems 
could be reflecting a systematic variation in gas-dynamical 
timescales as a function of galaxy mass and/or density. 

In any case, the fact that dynamical friction is unable to 
account for the steepening of an initially universal mass func- 
tion as the mass of the host galaxy decreases, combined with 
the possible existence in young, relatively unevolved cluster 
systems of a mass scale similar to Mc in our old GC systems, 
leads us to favor the view that a significant part of the ob- 
served morphology at the highest-mass ends of GCLFs is due 
to systematics in the initial distributions. The precise extent to 
which this part of the initial mass function is still reflected in 
the present-day one is still something of an open question, the 
answer to which will be a crucial ingredient in our understand- 
ing of GC formation and evolution. A detailed understanding 
of the "microscopic" star-formation processes on rather short 
timescales in very young clusters could well be key to making 
much further progress in this direction. 

8. SUMMARY AND CONCLUSIONS 

We have presented the GCLFs of 89 early-type galaxies 
in the Virgo cluster and determined maximum-likelihood es- 
timates for model parameters using fits of Gaussians and a 
simple "evolved Schechter function" described in 33.21 The 
latter reflects the effects of GC disruption (at a constant rate 
and presumably due mostly to two-body relaxation and evap- 
oration) on an initial cluster mass distribution that followed a 
Schechter function with a fixed power-law index of -2 at low 
masses. The evolved mass function tends to a flat shape at 
low M and is an accurate analytical approximation to the nu- 
merical distributions produced in the theory of Fall & Zhang 
(2001). We have tested the robustness of our results by sim- 
ulations, by the construction of GCLFs for galaxies binned 
together to contain a minimum number of clusters, and by us- 
ing alternate schemes to select GC candidates from catalogues 
of observed sources. Our main results and conclusions are the 
following: 

1. We find a remarkably regular decrease of the disper- 
sion of the GCLF as the luminosity of the host galaxy 
decreases (^and Jordan et al. 2006). Quantitatively, 
the maximum-likelihood estimates of the dispersion a 
of Gaussian fits to the z- and g-band data are well de- 
scribed by the linear relations presented in equations 
ST% and ( fTsT l. The dispersions for the GCLFs of the 
Milky Way and M3 1 fall in the midst of our new data 
and thus the correlation of cr with M^.gai would appear 
to be more fundamental than the older view, that GCLF 
widths depend on galaxy Hubble type. 

This trend reflects a systematic steepening of the GC 
mass function for massive clusters in particular (M > 
3 X 1O^M0, above the peak of the GCLF) as the host 
galaxy luminosity decreases. When fitting power-law 
mass functions to this upper cluster mass regime, the 
power-law exponents in a model of the form dN/dM cx 
increase from /3 < 2 to /? > 3 over the range of 
galaxy masses in our sample. This steepening is in turn 
equivalent to a systematic decrease of the cut-off mass 
Mc in evolved-Schechter function fits to the GCLFs, 
from Mc ~ 2-3 x 10^ Mq in the brightest galaxies to 
Mc ~ 3^ X 10^ M0 in the faintest systems. 

2. The GCLF turnover mass Mto is slightly smaller in 
dwarf systems (Mb -18), relative to the same quan- 
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tity in more massive galaxies. In the latter we have 
Mto = (2.2 ± 0.4) X lO^M©, decreasing to Mto ^ 1 -6- 
1 .7 X IO^Mq on average for the faintest galaxies in our 
sample — although individual dwarfs scatter between 
1 X 10^ Mq < Mto < 2 x 10^ Mq (®. We show that 
this might be at least partly accounted for by the effects 
of dynamical friction if all other processes shaping the 
mass function were to lead to an invariant Mto 07.21) . 

3. We explored radial variations of the GCLF over base- 
Unes of 20- 35 kpc in M87 ( VCC 1 3 1 6) and M49 ( VCC 
1226) by studying GCs in the fields of dwarf galaxies 
close in projection to these giant ellipticals ( ^6.3l l. We 
find no evidence for a variation of the turnover mass 
Mto with galactocentric distance in either galaxy, con- 
sistent with previous studies of M87 in particular. This 
reinforces the importance of the radial invariance of 
GCLFs as a constraint on models of GCLF formation 
and dynamical evolution. 

4. Our success in fitting evolved Schechter functions to 
our data 05.2I ) means that the GC mass functions in 
early-type Virgo galaxies are consistent with a univer- 
sally flat shape, dN/dM ^ constant, in the limit of low 
masses — as is also found in the Milky Way 03.3l and 
Fall & Zhang 2001). If this feature is caused by dy- 
namical evolution from a much steeper initial distri- 
bution, it requires that cluster masses decrease linearly 
in time. This can plausibly be expected if evaporation 
dominates the cluster evolution, although tidal shocks 
may also lead to similar behavior 

5. Fits of the evolved Schechter function imply that a nar- 
row range of average mass losses per GC — A « (2- 
10) X 1O^M0 at the outside — is required in all galax- 
ies to account for our observed GCLFs. Such a range 
of A across a factor of w 400 in galaxy luminosity is 
in rough agreement with observed (small) variations in 
the mean half-mass radii of GCs in the ACS VCS galax- 
ies (Jordan et al. 2005), and with simple scalings of 
evaporation rate as a function of host-galaxy luminos- 
ity 07.11 1. However, more work is required to reconcile 
fully the main idea — that long-term dynamical evolu- 
tion alone transformed initial Schechter cluster mass 
functions into the presently observed distributions — 
with the weak radial variation of GCLFs inside large 
galaxies and with observations of the orbital distribu- 
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tions and range of mean cluster densities in the same 
systems. 

6. The clear decrease of the GC cut-off mass with 
galaxy luminosity in evolved-Schechter function de- 
scriptions of the GCLF 05.2b is too pronounced to be 
explained by dynamical friction operating on a univer- 
sal dN/dM with an initially constant M^ in all galaxies 
07.2b . It most likely reflects systematic variations at 
the high-mass end of the initial GC mass function. 

The present-day mass functions of GCs were likely shaped 
by a variety of processes acting on different timescales, in- 
cluding systematic variations in the initial (proto-)cluster 
mass function at the high-mass end; long-term dynamical ero- 
sion by evaporation, tidal shocks, and dynamical friction; and 
global relaxation effects in time-varying galaxy potentials (hi- 
erarchical merging). It is further possible, though not yet en- 
tirely clear, that mass-selective early dissolution of clusters 
due to stellar evolution may have played some role in defin- 
ing the observed mass distributions. Future attempts to un- 
derstand the whole of the GCLF will clearly have to consider 
all of these processes, and their inevitable interplay, in quite 
some detail. Such comprehensive modeling will also have to 
acknowledge the increasingly complex and stringent empiri- 
cal constraints that follow from combining direct GCLF ob- 
servations with other GC systematics — such as their structural 
correlations, and the dynamics of cluster systems — for which 
data are continually accumulating and improving in quality. 
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0.275 


43 


1185 


15.68 


23.840 ±0.197 


0.691 ±0.137 


22.910 ±0.159 


0.639 ±0.1 13 


0.292 


33 




1489 


15.89 


23.977 ±0.439 


0.378 ±0.260 


23.156±0.381 


0.482 ±0.526 


0.417 


22 




1661 


15.97 


24. 177 ±0,154 


0.225 ±0.201 


23,059 ±0.417 


0,61 3 ±0,336 


0.477 


19 





Note. — Key to columns — (1) Galaxy VCC number; (2) Galaxy B magnitude from BinggeU et al. (1995); (3)-(4) Maximum-likelihood estimates of 
the Gaussian mean and dispersion a of the g-band GCLF; (5)-(6) Same as (3)-(4) but for the z-band; (7) Fraction B of the sample that is expected to be 
contamination; (8) Total number JV of all objects (including contaminants and uncorrected for incompleteness) with poc > 0-5; (9) Comments on individual 
galaxies. 



TABLE 3 

Definition of Binned GC Samples and Best-Fit GCLF Parameters 



Gaussian Fits 



Evolved Schechter Function Fits 



3roup 


N 1 


(Mo A 


B.gal 


B.gal 
















'-'z 


flic ' 






(j) 


(4) 


(j) 


(D) 


(') 


(°) 




(lU) 


(11) 


(12) 


U Jl 


(14) 





1 (1226) 


-21.8 


-21.8 


-21.8 


746 


-7.025 ± 0.086 


1.366 ±0.061 


-8.341 ±0.077 


1.321 ±0.053 


-7. 150 ±0.133 


-10.045 ±0.362 


-8.465 ±0.132 


-11.257 ±0.360 


1 


1 (1316) 


-21.5 


-21.5 


-21.5 


1721 


-7. 104 ±0.049 


1.312±0.035 


-8.433 ±0.041 


1.242 ±0.030 


-7.287 ±0.089 


-9.850 ±0.232 


-8.690 ±0.092 


-10.911 ±0.232 


2 


1 (1978) 


-21.3 


-21.3 


-21.3 


789 


-7.014 ±0.077 


1.340 ±0.058 


-8.329 ±0.070 


1.3 16 ±0.050 


-7.265 ±0.137 


-9.750 ±0.356 


-8.617 ±0.146 


-10.928 ±0.381 


3 


1 (881) 


-21.2 


-21.2 


-21.2 


355 


-7.334 ±0.097 


1.274 ±0.075 


-8.450 ±0.093 


1.238 ±0.071 


-7.533 ±0.198 


-9.877 ±0.525 


-8.607 ±0.221 


-11.043 ±0.647 


4 


1 (763) 


-21.1 


-21.1 


-21.1 


488 


-7.385 ±0.074 


1.178 ±0.055 


-8.522 ±0.070 


1.159 ±0.052 


-7.786 ±0.201 


-9.371 ±0.460 


-8.955 ±0.200 


-10.499 ±0.457 


5 


1 (731) 


-21.3 


-21.3 


-21.3 


888 


-7.431 ±0.061 


1.207 ±0.046 


-8.623 ±0.059 


1.199 ±0.044 


-7.651 ±0.137 


-9.789 ±0.339 


-8.818±0.134 


-11.011 ±0.337 


6 


1 (1903) 


-20.1 


-20.1 


-20.1 


294 


-7.434 ±0.089 


1.192 ±0.071 


-8.625 ±0.089 


1.215 ±0.073 


-7.641 ±0.209 


-9.971 ±0.611 


-8.776 ±0.188 


-11.418 ±0.576 


7 


1 (1632) 


-20.3 


-20.3 


-20.3 


439 


-7.089 ±0.103 


1.423 ±0.077 


-8.323 ±0.095 


1.390 ±0.071 


-7.198 ±0.152 


-10.393 ±0.474 


-8.443 ±0.154 


-11.516 ±0.469 


8 


1 (1231) 


-19.8 


-19.8 


-19.8 


239 


-7.221 ±0.090 


1.103±0.072 


-8.344 ±0.089 


1.106 ±0.069 


-7.841 ±0.311 


-8.888 ±0.669 


-9.013 ±0.319 


-10.002 ±0.680 


9 


2 


-19.6 


-19.7 


-19.5 


347 


-7.196 ±0.074 


1.102±0.057 


-8.308 ±0.076 


1.1 03 ±0.057 


-7.749 ± 0.279 


-8.909 ±0.635 


-8.946 ±0.294 


-9.957 ± 0.647 


10 


2 


-19.4 


-19.5 


-19.2 


248 


-7.282 ±0.088 


1.111 ±0.069 


-8.444 ± 0.089 


1.103 ±0.069 


-8.276 ±0.424 


-8.597 ±0.838 


-9.726 ±0.557 


-9.586 ± 1 .040 


11 


2 


-19.4 


-19.4 


-19.4 


283 


-7.341 ±0.086 


1.092 ±0.066 


-8.426 ±0.089 


1.096 ±0.068 


-8.485 ±0.517 


-8.478 ±0.974 


-9.431 ±0.453 


-9.696 ±0.879 


12 


2 


-19.1 


-19.3 


-18.9 


347 


-7.380 ±0.072 


1.091 ±0.055 


-8.447 ±0.071 


1.092 ±0.055 


-8.262 ±0.323 


-8.767 ±0.652 


-9.410 ±0.336 


-9.814 ±0.667 


13 


2 


-19.0 


-19.0 


-19.0 


212 


-7.333 ±0.084 


0.982 ± 0.065 


-8.446 ±0.083 


0.968 ± 0.066 


-8.321 ±0.519 


-8.467 ±1.013 


-9.596 ±0.602 


-9.482± 1.131 


14 


2 


-19.1 


-19.1 


-19.0 


214 


-7.459 ±0.084 


1.051 ±0.065 


-8.479 ±0.085 


1.042 ±0.068 


-8.028 ±0.335 


-9. 134 ±0.766 


-9.212 ±0.443 


-9.943 ± 0.964 


15 


4 


-18.5 


-18.6 


-18.3 


283 


-6.908 ± 0.088 


1.042 ±0.067 


-7.996 ±0.086 


1.032 ±0.066 


-8. 117 ±0.643 


-7.887 ±1.186 


-9.330 ±0.806 


-8.931 ±1.463 


16 


5 


-18.3 


-18.5 


-18.1 


257 


-7.081 ±0.089 


1.012 ±0.070 


-8. 146 ±0.087 


0.969 ± 0.068 


-8.399 ±0.648 


-7.972 ±1.158 


-9.413 ±0.642 


-9.057 ± 1.157 


17 


4 


-18.2 


-18.3 


-18.0 


208 


-7.338 ±0.087 


0.949 ±0.072 


-8.330 ±0.087 


0.945 ± 0.068 


-9.619±2.113 


-7.835 ±6.880 


-9.929 ±0.949 


-9.095 ±1.615 


18 


3 


-17.8 


-18.0 


-17.6 


205 


-7.263 ± 0.082 


0.951 ±0.062 


-8.276 ±0.085 


0.961 ±0.065 


-8.610 ±0.675 


-8.111 ±1.202 


-9.882 ±0.802 


-9.041 ±1.371 


19 


3 


-17.7 


-17.8 


-17.7 


197 


-7.409 ±0.079 


0.901 ±0.059 


-8.421 ±0.081 


0.905 ±0.060 


-10.101 ±1.986 


-7.839 ±7.800 


-10.042 ±0.896 


-9.133 ±1.512 


20 


8 


-17.1 


-17.5 


-16.8 


196 


-7. 149 ±0.099 


0.953 ±0.080 


-8.216 ±0.094 


0.927 ±0.072 


-8.247 ±0.680 


-8.025 ±1.253 


-9.371 ±0.630 


-9.161 ±1.154 


21 


6 


-16.6 


-16.8 


-16.5 


222 


-7.217 ±0.080 


0.943 ±0.060 


-8.240 ±0.080 


0.916 ±0.060 


-8.343 ±0.591 


-8. 166 ±1.103 


-9.609 ±0.656 


-9.079 ±1.155 


22 


9 


-16.4 


-16.7 


-16.1 


193 


-7.072 ±0.086 


0.875 ±0.068 


-8.043 ±0.096 


0.921 ±0.071 


-9.383 ±1.795 


-7.437 ±4.766 


-8.974 ±0.505 


-9.135 ±0.951 


23 


10 


-15.7 


-16.0 


-15.4 


201 


-7. 133 ±0.072 


0.749 ±0.058 


-8.090 ±0.077 


0.762 ±0.062 


-9.792 ±0.088 


-7.292 ±1.956 


-10.815 ±0.092 


-8.315 ±5.536 



o 

D 



Note. — Key to columns — (1) Identification number of binned group; (2) Number of galaxies that were used in the creation of this binned group. When only one galaxy present il 
st-fit Gaussian GCLF parameters and ct in the g- andz-bands; (11)-(14) Best-fit evolved Schechter GCLF parameters 5 and mc , in the g- andz-bands. 



VCC identifier is indicated; (3)-(5) Average, minimum, and maximum Afg of galaxies in this binned group; (6) Number of expected GCs; (7)-(10) 
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TABLE 4 
Best-fit power law exponent (3 



vcc 


0z 




VCC 


0z 




VCC 




Pg 


(1) 


(2) 


(3) 


(1) 


(2) 


(3) 


(1) 


(2) 


(3) 


1226 


1.80±0.11 


1.72±0.11 


654 


2.58 ±0.76 


2.73 ±0.72 


1178 


2.42 ± 0.40 


1.82 ±0.38 


1316 


1.75 ±0.07 


1.79 ±0.07 


944 


2.50 ±0.44 


2.55 ±0.46 


1283 


2.85 ±0.61 


2.90 ±0.61 


1978 


1.84±0.11 


1.77±0.11 


1938 


2.65 ±0.45 


2.83 ±0.45 


1261 


1.63 ±0.55 


2.01 ±0.52 


881 


1.80±0.16 


1.79±0.16 


1279 


1.87 ±0.28 


1.85 ±0.28 


698 


2.42 ±0.35 


2.38 ±0.33 


798 


2.15±0.15 


1.95 ±0.15 


1720 


2.69 ±0.47 


2.97 ±0.50 


1422 


2.95 ±0.94 


4.07±1.17 


763 


1.85±0.14 


1.87±0.14 


355 


3.50 ±0.97 


2.75 ±0.79 


2048 


1.44 ±0.87 


1.26±0.82 


731 


1.71 ±0.10 


1.77±0.10 


1619 


2.46 ±0.72 


2.25 ±0.71 


9 


2.42 ±0.82 


2.20 ±0.77 


1535 


2.03 ±0.20 


1.94 ±0.20 


1883 


3.18±0.60 


2.85 ±0.56 


1910 


2.35 ±0.52 


2.07 ±0.51 


1903 


1.87±0.18 


2.03 ±0.18 


1242 


3.25 ±0.45 


3.06 ±0.45 


856 


1.84 ±0.63 


1.70 ±0.64 


1632 


1.89±0.16 


1.83±0.16 


784 


3.77 ±1.14 


3.23 ±1.01 


140 


3.55 ±1.31 


2.51 ±1.21 


1231 


2.22 ±0.23 


2. 13 ±0.23 


1537 


2.13 ±0.62 


2.48 ±0.67 


1087 


2.73 ±0.54 


2.50 ±0.52 


2095 


1.79 ±0.34 


1.85±0.33 


778 


2.07 ±0.48 


2.00 ±0.48 


1861 


2.52 ±0.71 


1.92 ±0.62 


1154 


1.81 ±0.28 


2.02 ±0.28 


1321 


3.67 ±1.29 


4.99 ±1.97 


1431 


2.55 ±0.57 


2.69 ±0.59 


1062 


2. 13 ±0.26 


1.99 ±0.26 


828 


2.28 ±0.45 


2.48 ±0.42 


1528 


2.19 ±0.70 


2.56 ±0.67 


2092 


2.30 ±0.41 


2.35 ±0.41 


1250 


1.92 ±0.52 


2.09 ±0.49 


437 


3.55 ±0.91 


4.04 ±1.10 


369 


2.14±0.25 


2.13 ±0.25 


1630 


1.91 ±0.52 


1.99 ±0.50 


2019 


2.22 ±0.70 


3.17±0.81 


759 


2.32 ±0.27 


2.19 ±0.27 


1146 


2.33 ±0.52 


2.47 ±0.50 


21 


2.32 ±0.88 


2.76 ±0.93 


1692 


1.93 ±0.29 


2.40 ±0.29 


1025 


3.09 ±0.75 


2.81 ±0.66 


1499 


2.75 ±0.90 


2.69 ±0.88 


1030 


1.93 ±0.25 


2.13 ±0.25 


1303 


2.55 ±0.59 


2.48 ±0.54 


1545 


2.61 ±0.74 


2.57 ±0.73 


2000 


2.07 ±0.25 


2. 16 ±0.24 


1913 


3.03 ±0.58 


2.57 ±0.58 


1075 


3.94± 1.31 


3.85 ± 1.28 


685 


1.71 ±0.25 


1.71 ±0.24 


1125 


2.78 ±0.57 


2.37 ±0.58 


1539 


3.14±0.92 


2.69 ±0.85 


1664 


1.85 ±0.29 


1.66 ±0.28 


1475 


2.37 ±0.54 


2.55 ±0.55 


1185 


5.63 ±1.62 


5.56±1.59 



Note. — Key to columns — (1) Galaxy VCC number; (2) Best-fit power-law exponent {3 for the mass function of GCs between ~0.5-2.5 
mag brighter than the turnover magnitude in the z-band GCLF; (3) As (2) but for the g-band. 
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TABLE 5 

Average GC Colors and Mass-to-Light Ratios 



vcc 




T- 




VCC 


/(0 — 7'\\ 

\ys <->) 




T 


VCC 




T- 


T 

^g 








y^f 


y^j 


y^f 


y-') 


(A) 

y^) 


y^) 


y^) 




(41 

y^} 


1226 


1.24 


1.47 


2.69 


1242 


1.11 


1.49 


2.44 


1297 


1.05 


1.50 


2.33 


1316 


1.23 


1.47 


2.67 


784 


1.14 


1.48 


2.50 


1861 


1.00 


1.50 


2.24 


1978 


1.25 


1.47 


2.72 


1537 


1.00 


1.50 


2.24 


543 


0.94 


1.51 


2.12 


881 


1.09 


1.49 


2.41 


778 


1.04 


1.50 


2.31 


1431 


1.00 


1.50 


2.24 


798 


1.14 


1.48 


2.50 


1321 


1.04 


1.50 


2.31 


1528 


0.95 


1.51 


2.14 


763 


1.11 


1.49 


2.44 


828 


1.00 


1.50 


2.24 


1695 


1.01 


1.50 


2.26 


731 


1.19 


1.47 


2.59 


1250 


0.98 


1.51 


2.20 


1833 


1.01 


1.50 


2.26 


1535 


1.18 


1.48 


2.57 


1630 


1.10 


1.49 


2.42 


437 


0.90 


1.52 


2.05 


1903 


1.18 


1.48 


2.57 


1146 


1.20 


1.47 


2.61 


2019 


0.90 


1.52 


2.05 


1632 


1.21 


1.47 


2.63 


1025 


0.97 


1.51 


2.18 


200 


0.82 


1.54 


1.91 


1231 


1.12 


1.48 


2.46 


1303 


0.94 


1.51 


2.12 


571 


0.92 


1.52 


2.09 


2095 


1.07 


1.49 


2.37 


1913 


1.02 


1.50 


2.27 


21 


0.88 


1.52 


2.01 


1154 


1.12 


1.48 


2.46 


1327 


1.06 


1.49 


2.35 


1488 


0.87 


1.52 


1.99 


1062 


1.14 


1.48 


2.50 


1125 


0.93 


1.51 


2.11 


1499 


0.93 


1.51 


2.11 


2092 


1.13 


1.48 


2.48 


1475 


0.94 


1.51 


2.12 


1545 


0.93 


1.51 


2.11 


369 


1.15 


1.48 


2.52 


1178 


1.06 


1.49 


2.35 


1192 


1.10 


1.49 


2.42 


759 


1.10 


1.49 


2.42 


1283 


1.03 


1.50 


2.29 


1075 


0.93 


1.51 


2.11 


1692 


1.08 


1.49 


2.39 


1261 


1.05 


1.50 


2.33 


1440 


0.98 


1.51 


2.20 


1030 


1.14 


1.48 


2.50 


698 


1.00 


1.50 


2.24 


230 


0.92 


1.52 


2.09 


2000 


1.05 


1.50 


2.33 


1422 


1.09 


1.49 


2.41 


2050 


0.89 


1.52 


2.03 


685 


1.07 


1.49 


2.37 


2048 


1.01 


1.50 


2.26 


751 


0.85 


1.53 


1.96 


1664 


1.18 


1.48 


2.57 


1871 


0.96 


1.51 


2.16 


1828 


0.88 


1.52 


2.01 


654 


0.99 


1.50 


2.22 


9 


1.01 


1.50 


2.26 


1407 


1.02 


1.50 


2.27 


944 


1.06 


1.49 


2.35 


575 


1.00 


1.50 


2.24 


1886 


0.80 


1.55 


1.90 


1938 


0.99 


1.50 


2.22 


1910 


1.06 


1.49 


2.35 


1199 


1.13 


1.48 


2.48 


1279 


1.04 


1.50 


2.31 


1049 


0.97 


1.51 


2.18 


1539 


0.97 


1.51 


2.18 


1720 


1.08 


1.49 


2.39 


856 


1.02 


1.50 


2.27 


1185 


0.92 


1.52 


2.09 


355 


1.09 


1.49 


2.41 


140 


1.00 


1.50 


2.24 


1489 


0.98 


1.51 


2.20 


1619 


1.06 


1.49 


2.35 


1355 


0.92 


1.52 


2.09 


1661 


0.95 


1.51 


2.14 


1883 


1.06 


1.49 


2.35 


1087 


0.94 


1.51 


2.12 











Note. — Key to columns — (1) Galaxy VCC number; (2) M ean GC (g-z) color (from data in Peng et al. 2006a); (3)-(4) Average GC 
mass-to-light ratio in the z and g bands, obtained as described in § 16.1.21 



TABLE 6 

Gaussian GCLF Parameters for Outer GCs of M87/M49 Companions' 



VCC 




"g 






B 


N 


Host 




Tg 




(I) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


1327 


23.886±0.182 


1.288 ±0.144 


22.777 ±0.164 


1.224±0.129 


0.070 


93 


VCC1316 


1.104 


2.43 


1.49 


1199 


23.805±0.127 


1.224 ±0.102 


22.631 ±0.117 


1.175±0.091 


0.053 


151 


VCC 1226 


1.100 


2.42 


1.49 


1192 


23.643 ±0.102 


1.013 ±0.080 


22.540 ±0.099 


1.011 ±0.076 


0.064 


144 


VCC 1226 


1.155 


2.54 


1.48 


1297 


23.410±0.183 


1.208 ±0.142 


22.320 ±0.166 


1.138±0.130 


0.103 


67 


VCC1316 


1.090 


2.40 


1.49 



Note. — Key to columns — (1) Galaxy VCC number; (2)-(3) Maximum-likelihood estimates of the Gaussian mean ^ and dispersion cr 
of the g-band GCLF; (4)-(5) Same as (2)-(3) but for the z-band; (6) Fraction B of the sample that is expected to be contamination; (7) Total 
number A' of all objects with pac > 0.5 (including contaminants and uncorrected for incompleteness); (8) VCC number of giant elliptical 
galaxy close in projection; (9) Mean GC ig-z) color; (lO)-(l 1) Average GC mass-to-light ratio in the g and z bands, obtained as described in 
§11X21 

" All reported numbers refer only to those GC candidates that are more than 6 {Re) away from the centers of the galaxies indicated, where 
{Rc) is the median effective radius of other VCS galaxies that have magnitudes within 0.5 mag of the target galaxy. 



